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CHAPTER  I 
INTRODUCTION 


With  the  advent  of  complex  equipment  in  both  the  government 
and  industrial  sectors  of  the  economy,  the  study  of  system 
reliability  and  maintainability  has  become  increasingly  important. 
Complex  systems,  such  as  aircraft  and  mortar  artillery  locating 
radar  systems,  must  frequently  operate  over  long  periods  of  time 
far  from  their  sources  of  supply.  It  often  happens  that  spares 
of  a given  part  type  are  exhausted,  while  spares  of  another  part 
type  are  in  excess;  subsequent  failures  of  the  part  type  with 
a deplete  spare  inventory  will  terminate  system  operation  until 
replacements  arrive  from  some  external  source  of  supply.  This 
anomoly  is  a source  of  exhorbitant  logistics  support  costs,  when 
the  cost  of  lost  service  or  "downtime"  is  considered  in  the 
overall  cost  model  of  the  system,  for  a system  operating  in  the 
private  industrial  sector,  the  loss  of  system  performance  may 
alternatively  be  penalized  with  the  dollar  cost  of  lost  revenue. 

The  potential  cost  of  such  lost  service  provides  the  motivation 
for  developing  a working  policy  for  maintaining  a particular  system. 
The  specific  technique  studied  in  tin's  research  is  cannibalization. 

The  citations  on  the  following  pages  follow  the  style  of  the 
Naval  Re  sea  roll  Logistics  Qua  rterl  y_. 
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Cannibalization,  as  defined  in  this  context,  is  the  replenish- 
ment of  the  spare  inventory  with  working  components  from  a 
terminally  failed  unit. 

For  purposes  of  illustration,  consider  a system  made  up  of 
N independent  and  identical  units,  each  composed  of  P serial 
components.  When  a component  fails,  it  is  replaced  with  a spare 
component  from  the  inventory;  if  there  are  no  available  spares  for 
this  particular  component  type  the  unit  is  disassembled,  the 
remaining  working  components  are  restored  to  original  condition 
and  added  to  their  respective  spare  inventory  locations.  Thus, 
cannibalization  is  the  procedure  of  disassembling  the  remaining 
working  components,  restoring,  and  adding  them  to  the  remaining 
spare  component  types.  In  effect,  it  is  the  dismantlement  of 
units  for  components  to  be  used  as  replacements  in  other  units. 

It  must  be  realized,  though,  that  this  technique  is  appropriate 
and  useful  only  when  specific  logistics  problems,  such  as  very 
large  procurement  lead  times,  prevent  the  acquisition  of  a 
replaceable  component  from  the  primary,  or  an  alternate,  source  of 
supply.  Otherwise  it  would  obviously  be  appropriate  to  wait  for  the 
arrival  of  the  the  failed  component  type. 

The  traditional  concept  of  cannibalization  is  one  in  which  the 
working  components  of  a failed  unit  are  utilized  directly  as 
replacements  for  other  units  as  their  components  successively  fail. 

The  primary  difference  between  the  "traditional"  definition  of 
cannibalization,  and  the  one  used  here,  is  that  in  the  traditional 


! 
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definition  the  replacement  activity  must  encompass  both  the  dis- 
assembly from  the  donor  system  and  subsequent  installation  into 
the  failed  system;  whereas  in  the  more  restrictive  definition 
used  in  this  application,  the  replacement  activity  is  defined 
strictly  as  a spare  extraction  from  an  inventory  of  discrete 
component  types.  The  cannibal ization  activity  is  quantitatively 
defined  in  such  a way  to  adequately  describe  the  stochastic 
process  of  terminal  component  failure. 

The  tradition!  concept  is  not  addressed  in  tin's  research 
because  it  requires  a more  complex  definition  of  the  replacement 
time  density  function.  Specifically,  the  replacement  time  of  a 
terminally  failed  component  within  a working  system  would  have  to 
include  the  time  required  to  locate  and  remove  the  correct  part 
from  the  donor  system  as  well  as  the  time  required  for  installation 
into  the  recipient  system.  The  more  restrictive  definition  of 
cannibalization  gives  a simpler  stochastic  definition  of  the 
probability  of  cannibalizing  any  unit  because  of  the  failure  of  a 
particular  component  type. 

The  technique  of  cannibalization  can  provide  savings  in  the 
operational  costs  of  a system  because  it  limits  or  reduces  the 
required  number  of  initial  spares  for  a working  system;  also, 
cannibalization  offers  a method  for  "re-using"  most  of  the 
components  of  a terminally  failed  unit,  which  otherwise  would  have 
a small  or  negligible  salvage  value. 

In  today's  competitive  economic  situation,  doubly  plagued 


by  the  simultaneous  effects  of  inflation  and  recession,  it  is  not 
economically  feasible  to  overdesign  with  respect  to  system  reliabil 
ity,  maintainability,  or  availability.  An  optimal  economic  system 
design  should  yield  a minimal  cost  while  still  meeting  a minimum 
availability  or  reliability  requirement.  The  system  availability, 
A(t),  is  defined  as  the  probability  that  a system  is  functioning 
at  any  given  point  in  time,  t.  On  the  other  hand,  the  system 
reliability,  R(t),  is  the  probability  that  a system  has  functioned 
successfully,  that  is,  it  has  experienced  no  failures  for  an  in- 
terval of  time,  [0,t]  . Thus,  in  general,  R(t) 

Another  feasible  model  constraint  is  a requirement  that  the 
probability  of  experiencing  a critical  spare  shortage,  which 
could  require  cannibalization  of  a particular  unit,  be  strictly 
less  than  some  preassigned  value,  say  1-a.  This  constraint  can  be 
stated  alternatively  as  follows:  the  probability  of  no  critical 
spare  shortage  for  the  system  must  be  at  least  as  large  as  ex.  By 
"critical  spare  shortage"  is  meant  an  inventory  depletion  which 
precludes  further  system  operation  upon  subsequent  failure  of  any 
component  type. 

It  is  the  purpose  of  this  research  to  quantitatively  analyze 
cannibalization  as  a maintainability  technique  for  minimizing  the 
present  worth  of  system  operational  costs  over  the  system  lifetime 
subject  to  some  constraint  on  system  performance. 

It  is  emphasized  that  the  selection  and  application  of  the 
proper  maintainability  technique  results  iri  nwny  economies  measured 


in  terms  of  labor  time,  material,  and  money.  These  savings  are 
attributable  to  the  fact  that  maintainability  prediction  is  conside 
ed  to  be  a tool  for  design  enhancement  because  it  provides  for  the 
early  recognition  and  elimination  of  areas  of  poor  maintainability 
during  the  early  stages  of  the  design  life  cycle.  Otherwise,  areas 
of  poor  maintainability  would  only  become  apparent  during  demon- 
stration testing  or  actual  use,  after  which  time  correction  of 
design  deficiencies  would  be  costly,  and  unduly  delay  schedules 
and  missions. 

Besides  simply  being  a maintainability  technique,  this 
formulation  of  cannibalization  may  be  implemented  as  a design  tool 
for  the  initial  system  spare  configuration  which  would  minimize 
the  overall  operating  costs  to  the  user.  It  is  especially 
appropriate  to  systems  which  experience  deterioration  and  ob- 
solescence with  time  as  was  evidenced  with  the  B- 52  and  DC-3 
aircraft.  A functional  policy  of  cannibalization  may  have  been 
able  to  prolong  the  useful  life  of  such  aircraft  in  Southeast 
Asia  with  a smaller  expense  to  the  Defense  Department. 
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Research  Objectives  and  Assumptions 

The  primary  goal  of  this  research  is  to  develop  an  operational 
cost  model  which  allows  one  to  consider  cannibalization  as  a main- 
tainability technique.  The  analysis,  based  on  specific  component 
failure  and  replacement  density  functions,  is  a systematic  presen- 
tation of  the  cannibalization  technique,  and  its  consequences  in 
determining  the  optimal  spare  inventory  allocation  for  a system  of 
N independent  and  identical  units,  each  with  P discrete  components 
or  modules  per  unit.  The  optimal  spare  inventory  allocation  is  the 
one  with  the  minimal  present  worth  of  the  net  operational  costs,  yet 
satisfies  a specific  constraint  on  overall  performance,  such  as 
meeting  a minimum  level  of  overall  availability,  or  a- requirement 
that  the  probability  of  no  critical  spare  shortage  for  the  system 
be  at  least  as  large  as  a pre-assigned  value,  say  ao,  specified 
by  the  system  user.  Because  the  number  of  spares  are  finite, 
a constraint  on  system  availability  is  not  mathematically  feasible, 

except  in  the  steady  state  when  the  availability  reaches  zero  for 

0 

a finite  number  of  units  and  spares. 

A possible  mathematical  formulation  of  the  problem  is  stated 
as  follows: 

T N 

Minimize  K = K + K?  + K,  - IT.  V At)  P.  (t)  exp(-rt)  dt, 

0 1 c 6 0 i=0 

where  Kq  is  the  present  worth  of  the  system  operational  costs,  K,  is 
the  present  worth  of  the  manufacturing  or  procurement  costs,  is 


the  present  worth  of  the  replacement  costs,  and  is  the  present 
worth  of  the  cannibalization  costs.  V.(t)  is  the  system  return  or 
revenue  function  when  i units  are  operating  at  time  t,  P^t)  is  the 
probability  of  exactly  i units  operating  at  time  t. 

The  present  worth  of  the  return  or  revenue  function  is  given  by 

the  integral  of  the  product  of  the  expected  return  for  having  at 
least  0 units  operating,  with  the  continuous  compounding  discount 
factor,  exp(-rt).  The  present  worth  calculation  is  performed  over 
the  closed  interval  [o,f|  , where  T is  the  specified  system  time 
horizon. 

The  development  of  the  model  will  be  subject  to  the  following 
assumptions: 

1.  Each  of  the  f!  units  is  identical,  and  consists  of  P 
components. 

2.  The  system  is  supported  with  an  inventory  of  spares,  with 

m.  spares  initially  stocked  for  the  j^  component  type, 

J 

1 ^ j j<  P.  Spares  are  not  subject  to  failure,. 

3.  All  failures  are  independent  and  the  time  to  failure  for 
the  j^1 2 3 4  component  type  is  governed  by  the  exponential 
probability  density  function,  f.(t)  = exp  (-X.jt),  for 
1 < j < P. 

4.  Failures  are  detected  instantly  and  replacement  commences 
immediately  if  a spare  for  the  particular  component  type 
is  available,  otherwise  the  unit  is  cannibalized. 


10.  System  operational  costs  consist  only  of  manufacturing  or 
procurement  costs,  replacement  costs,  cannibalization 
costs,  and  a revenue  or  return  function  which  is  treated 
as  a negative  cost. 

Dissertation  Outline 

Chapter  II  contains  a thorough  review  of  the  currently 
published  literature  on  cost  modelling  and  maintainability  analysis. 


Chapter  III  develops  the  overall  system  cost  model  for  a 
general  system  configuration  and  contains  a stochastic  description 
of  the  process  of  cannibalization.  It  results  in  an  explicit 

equation  for  the  cannibalization  probability  as  a function  of  time. 
Also  included  in  Chapter  III  are  some  example  cost  calculations 
and  a discussion  of  possible  revenue  or  return  functions. 

Chapter  IV  outlines  applicable  system  constraints  and 
illustrates,  in  detail,  the  required  calculations  for  a constraint 
on  spare  adequacy.  Other  constraints  discussed  are  limitations 
on  the  total  weight  or  volume  of  a specific  spare  allocation, 
a resource  constraint  on  the  total  investment  in  spares,  a 
requirement  on  a system's  interval  availability,  and  a maintain- 
ability constraint  on  system  downtime. 

Chapter  V contains  the  necessary  criteria  for  applying  this 
cost  model  and  analyzes  an  example  problem  in  detail.  Also 
included  in  Chapter  V is  the  State  Transition  Analysis  necessary 
to  calculate  P-j ( t ) , the  probability  of  exactly  i units  operating 
at  time  t. 

Chapter  VI  contains  a concluding  discussion,  and  suggestions 


for  further  research. 
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on  the  optimal  transition  of  multi -unit-component  systems  was  done 
by  Brown [7]  . The  optimal  transition  policy  was  the  one  which 
yielded  the  maximum  net  revenue  to  the  system  user  over  a prescribed 
system  lifetime.  The  policy  determined  if  and  when  to  retire  a 
system  from  active  use,  and  use  it  as  a source  of  spares  for  other 
fielded  systems.  In  this  work,  lie  developed  an  algorithm  based  on 
the  present  worth  of  the  system  return  as  the  decision  criterion. 

He  did  not  consider  the  overall  system  availability  or  the 
probability  of  spare  adequacy  as  constraints  in  his  model. 

Brown  made  nine  general  assumptions  in  his  model: 

1.  The  spares  are  not  subject  to  failure. 
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2.  All  the  components  of  a unit  must  be  functioning  for  a 
unit  to  be  available. 

3.  The  terminal  failures  of  the  components  are  statistically 
independent. 

4.  The  replacement  or  cannibalization  of  a component  does 
not  cause  failures. 

5.  Replacement  times  for  failed  components  are  sufficiently 
short  to  not  affect  the  availability  of  a unit.  (This  is 
equivalent  to  an  instantaneous  replacement  model.) 

6.  The  parameters  required  by  the  model  can  be  determined. 

7.  Only  random  failures  of  the  components  are  considered. 

8.  Failures  are  amenable  to  the  Poisson  postulates  thus 
having  an  exponential  distribution  of  time  to  failure. 

9.  Salvage  values  are  negligible. 


Although  Brown  ackrtowl edged  that  there  existed  a time  depend- 
ent probability  of  cannibalization,  as  well  as  a time  dependent 
probability  of  replacement,  he  never  developed  a closed  form 
description  of  these  functions. 

Brown  also  considered  the  incremental  changes  in  the  system 
return  function  due  to  incremental  changes  in  the  number  of  spares 
for  any  component  type.  He  also  discussed  the  effects  of  perturbing 
the  number  of  units  on  the  system  return  function.  Both  analyses 
were  done  using  mathematical  induction  rather  than  any  generalized 
closed  form  analysis  or  system  simulation. 

Gaver  [l 5]  studied  the  availability  of  some  simple  redundant 
systems,  in  which  repair  was  possible  in  case  of  failure.  The  systems 
discussed  were  composed  of  two  identical  subsystems,  and  the  system 
was  considered  to  be  in  a state  of  failure  when  both  subsystems  were 
simultaneously  in  such  a state.  Gaver1 s work  was  significant  in  that 
it  showed  the  relative  gains  possible  in  the  mean  time  to  system 


failure  when  repair  facilities  are  added, 
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Wi 1 ken  and  Langford  [51]  determined  the  probability  of  exactly 
i failures  occurring  in  a time  period  of  length  t,  P(i,m,N,t),  for  a 
system  of  N identical  single  component  units  in  parallel  with  m 
spares.  They  derived  the  following  equation,  which  is  correct  for 
exponential  failures  and  instantaneous  replacement. 

f(NAt)^  exp(-NAt)/i ! , for  0 < i < m. 


P(i,m,N,t)  =< 


Nm  N! 


]i-m) ! (N-i+m) , 


i-m-1 

Z 

0=0 
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1)J  (i-m) 
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exp( (i-m-j)At)  _ 
(i-m-j)m 


+m) ! 
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Z (i-m-j) 
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for  i+m  < i < m+N. 
for  i > m+N 


Guillory  [18]  illustrated  a policy  of  backrobbing,  in  which  a 
unit  or  sub-system  is  removed  from  the  repair  queue  in  order  to 
expedite  the  repair  of  a failed  system  inside  the  service  facility. 
The  technique  was  used  when  the  inventory  of  spares  was  deficient 
and  the  repair  time  would  cause  delay  in  the  repair  of  the  system. 

Guillory's  analysis  was  important  primarily  from  the  standpoint 
of  expediting  the  repair  of  critical  systems.  It  emphasized  the 
stochastic  behavior  of  systems  within  a repair  queue;  however  it  did 
not  fully  address  the  aspect  of  operational  costs  and  spares 
provisioning,  as  is  intended  by  this  research. 


Hirsch,  Meisner  and  Boll  [19]  developed  a "structure  function" 
in  order  to  obtain  an  effective  method  for  determing  the 
probability  density  function  of  the  performance  levels  of  a system 
subject  to  cannibalization.  They  assumed  that  parts  either  "work" 
or  "fail",  but  allowed  the  possibility  that  the  system  may  function 
at  any  one  of  several  levels  of  performance. 

They  studied  the  probability  distribution  of  system  states  in 
the  two  special  cases: 

1.  No  spares  are  available. 

2.  Spare  part  inventories  are  arbitrary,  and  failures  of 
parts  are  governed  by  the  Poisson  Process. 

They  also  made  the  following  two  assumptions: 

1.  Failures  are  detected  instantaneously,  and  part  replace- 
ments and  interchanges  a>~e  performed  instantaneously. 

2.  Part  lifetimes  are  independent  random  variables,  and  parts 
of  the  same  type  have  identical  lifetime  distributions 
which  are  not  affected  by  cannibalization. 

Again,  the  work  of  Hirsch,  Meisner,  and  Boll  was  primarily  one 
of  developing  a particular  stochastic  description  of  the  process  of 
cannibalization  without  regard  to  the  potential  savings,  in  terms  of 
decreased  system  operational  costs  throughout  the  life  of  the  system. 

Simon  [40,  4l]  generalized  the  model  of  Hirsch,  Meisner  and  Boll 
and  developed  an  "admissible"  cannibalization  policy  which  allocated 
shortages  to  maximize  the  system  state  following  each  failure.  He 
showed  that  not  all  admissible  policies  are  equally  desirable. 

Simon's  work  was  essentially  an  extension  of  the  model 
described  by  Hirsch,  Meisner,  and  Boll,  with  no  emphasis  on  the 
aspects  of  overall  system  costs  or  the  provisioning  of  spares. 


Similarly,  the  research  results  presented  in  this  Dissertation 
are  based  upon  the  present  worth  of  particular  operational  costs  as 
the  decision  criterion  for  the  system  configuration.  However,  no 
generalized  operational  cost  patterns  are  assumed;  rather,  an 
expression  for  the  present  worth  of  the  system  operational  costs  is 
derived.  In  the  analysis  presented  herein,  the  assumed  revenue 
function  is  treated  as  a negative  cost. 

McNichols  and  Messer  [23]  developed  a procedure  for  allocating 
"availability  parameters  (repair  times  and  failure  rates)"  to  the 
specific  components  of  series  systems  in  order  to  minimize  the  cost 
of  system  development.  They  used  the  Lagrange  Multiplier  technique 
of  constrained  optimization  assuming  constant  failure  rates  and 
repair  times.  Their  cost  function  specifically  delineated  the  cost 
incurred  in  order  to  decrease  the  failure  rates  of  each  of  the 
components  as  well  as  the  cost  of  decreasing  the  repair  times.  The 
authors  labeled  these  costs  "component  improvement  costs". 

They  solved  for  the  optimal  system  parameters  of  failure  rate 
and  repair  time  for  eacii  of  the  respective  components  which 
minimized  the  system's  improvement  cost  required  to  meet  a particular 
availability  level.  They  also  investigated  the  effects  of  holding  a 
parameter  (failure  rate)  at  a fixed  value  while  allowing  the 
parameters  of  the  remaining  components  to  change.  They  pointed 
out  that  the  constant  failure  rate  assumption  was  necessary  in 
order  to  develop  and  optimize  the  analytic  model. 

They  did  not  investigate,  however,  the  effects  of  sparing  and 
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downtime  upon  the  overall  system's  cost,  because  the  model  was  to  be 
implemented  during  on  equipment's  design  stage  when  knowledge  of 
sparing  and  downtime  costs  is  limited. 

Morrison  [29]  considered  a complex  system  whose  components  were 

divided  into  two  subsystems,  each  containing  components  whose  lives 

were  exponentially  distributed  with  different  scale  parameters.  The 

components  were  assumed  connected  in  series,  thus  failure  of  any 

component  caused  the  entire  system  to  fail.  Each  of  the  subsystems 

was  unique  and  supplied  with  an  inventory  of  spares.  System 

fail  ures  were  corrected  by  replacing  failed  subsystem  components 

with  components  from  the  appropriate  spare  inventory.  When  a 

; 

failure  occurred  with  the  corresponding  spare  inventory  deplete,  the 
entire  system  was  discarded.  No  form  of  cannibalization  was  con- 
sidered because  there  were  only  two  unique  subsystems,  each  with 
different  failure  density  parameters,  and  presumably,  different 
functional  characteristics.  He  gave  an  explicit  solution  to 
the  specific  problem  of  allocating  spares  to  maximize  the 
expected  life  of  a system  with  an  inventory  of  similar  spares. 

His  analysis  was  restricted  to  components  with  an  exponential 
failure  density.  He  developed  tables  of  the  maximum  ex- 
pected system  life  under  various  spare  allocations  for  each 
of  the  component  types.  He  also  allocated  spares  for  maximum 
rel iabi 1 ity. 

Again,  this  work  was  primarily  directed  to  allocating  spares  to 
maximum  reliability  without  considering  system  costs.  It  also 


neglects  the  density  function  of  time  to  repair. 

Morrison  and  David  [30j  obtained  the  distribution  of  operating 
life  for  a series  system  of  like  elements  supplied  with  a set  of 
spare  components.  This  distribution  has  been  evaluated  for  some 
common  types  of  component  life  densities,  exponential  and  gamma,  and 
tables  of  expected  total  system  life  have  been  constructed.  These 
exnectations  were  compared  with  those  of  systems  with  no  spares  as 
a measure  of  the  efficacy  of  the  additional  spare  components.  They 
also  studied  the  reliability  of  systems  with  spares. 

The  report  by  Morrison  and  David  is  similar  to  the  research  of 
this  Dissertation,  in  that  several  failure  distributions  were 
considered,  and  the  system  configuration  was  initially  a series 
system  of  identical  elements,  supplied  with  a set  of  spare 
components.  The  system  efficacy  was  measured  by  comparing  the  ex- 
pected system  life  for  a system  with  a particular  number  of  spares  to 
the  expected  life  of  a system  with  no  spares.  However,  as  in  most  of 
the  research  cited  so  far,  no  emphasis  was  placed  on  the  cost  of 
deploying  and  operating  a particular  system. 

Crow  [12]  considered  the  theoretical  and  practical  implications 
of  the  nonhomogeneous  Poisson  Process  model  for  reliability,  and 
gave  estimation,  hypothesis  testing,  comparison  and  goodness  0^  fit 
procedures  when  the  process  had  a Ueibull  probability  density 
function.  Applications  of  the  Weibull  model  in  the  field  of 
reliability  and  in  other  areas  were  discussed. 

Crow  postulated  that,  for  a repairable  system,  one  is  rarely 
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interested  primarily  in  the  time  to  first  failure.  Rather,  he 
emphasized  the  interest  in  the  probability  of  system  failure  as  a 
function  of  system  age.  He  pointed  out  that  often  mathematical 
idealizations  assume  that  the  failure  times  of  the  complex  repairable 
system  follow  a Poisson  Process. 

For  various  repairable  systems,  particularly  of  the  complex 
electronic  type,  a constant  failure  rate  had  been  considered 
representative,  after  perhaps  an  initial  burn  in  period.  The  period 
in  which  a system  exhibits  a constant  intensity  of  failure  is  often 
called  the  system  "useful  life".  If  a complex  system  consists  of  a 
large  number  of  components,  each  acting  independently,  if  the  failure 
of  a component  results  in  a failure  of  the  system,  and  if  each 
component  is  replaced  upon  failure,  then  under  fairly  general 
conditions,  the  occurrence  of  failures  of  the  system  will  approach  a 
constant  intensity,  or  rate,  as  the  number  of  components  and 
operating  time  become  large.  However,  many  large  or  complex 
repairable  systems  generally  experience  a wearout  phase  which 
eventually  makes  them  economically  impractical  or  too  unreliable  to 
continue  in  service  without  undergoing  overhaul.  Such  systems, 
according  to  Crow,  will  never  achieve  the  equilibrium  state  of  a 
homogeneous  Poisson  Process. 

Crow's  model  was  primarily  directed  to  a stochastic  system 
description,  rather  than  to  a cost  modelling  methodology  as  is  a 
primary  objective  of  this  research. 


Timsans  et  al . [45]  developed  a family  of  hyperbolic  cost 
functions  which  could  be  used  to  analyze  the  cost  of  achieving 
prescribed  levels  of  reliability  and  maintainability,  for  series 
systems  experiencing  constant  failure  rates.  The  required  failure 
rates,  and  mean  repair  times,  were  determined  in  order  to  minimize 
the  total  system  cost  yet  meet  an  availability  constraint.  They 
used  a Lagrange  Multiplier  technique  essentially  equivalent  to 
that  used  by  McNichols  and  Messer  [23]  to  solve  the  allocation 
problem. 

Cleroux  and  Hans com  [9]  analyzed  age  replacement  policies 
where  replacement  of  a unit  occur-  at  failure,  or  at  age  T, 
whichever  occurs  first.  The  age  replacement  policy  which  min- 
imized the  average  expected  cost  per  unit  time  over  an  infinite 
time  span  was  obtained  in  t^  whore  the  cost  structure  in- 
volved a term  accounting  for  adjustment  costs  or  interest  charges 
which  were  incurred  at  fixed  intervals  of  time  of  equal  length. 

The  adjustment  costs  included  maintenance  expenses  required  to 
keep  the  equipment  operating,  and  depreciation  costs  resulting 
from  the  decline  in  resale  value  due  to  age  and  condition.  However 
no  consideration  was  given  to  the  density  of  time  to  perform  repair 
and  replacements. 

Proschan  [33]  solved  the  inventory  problem  of  maximizing  spare 
adequacy  subject  to  a single  linear  constraint  on  the  total  in- 
vestment in  spares.  He  modelled  the  process  with  exponential  fail- 
ures and  instantaneous  replenishment. 
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He  also  gave  a thorough  discussion  of  the  use  of  Polya  type 
distributions  in  Renewal  Theory.  He  proved  that  the  Polya  distri- 
bution family  is  closed,  i.e.  the  convolution  of  a Polya  type  m 
function  with  a Polya  type  n function  is  a Polya  function  of 
type:  Min(m,n). 

None  of  the  previous  works,  which  dealt  with  the  stochastic 
analysis  of  systems  subject  to  failure  and  repair,  included  the 
repair/replacement  density  function  in  the  analysis.  Instantaneous 
replenishment,  or  fixed  repair  times,  were  generally  assumed  in 
most  of  the  stochastic  models  developed.  Consequently,  a cost 
model  formulated  using  the  relevant  repair/replacement  density 
functions  is  a significant  contribution  to  the  state  of  the  art, 
and  extends  the  scope  of  the  current  literature.  The  analysis 
contained  herein  is  such  a contribution. 


The  author  develops  a stochastic  analysis  of  the  repair/replace- 
ment and  cannibalization  processes,  and  derives  general  formulas 
for  the  time  dependent  repair/replacement,  and  cannibalization 
probability;  the  expressions  are  based  upon  general  failure  and 
replacement  density  functions.  He  also  performs  a State  Transition 
Analysis,  using  an  Absorbing  State  Markov  Process,  to  model  the 
consumption  of  spare  components. 

The  cost  model  may  be  subject  to  one  or  more  constraints, 
which  provide  bounds  on  the  required  number  of  spares  for  each 
component  type. 

An  example  is  presented,  and  numerical  solutions  are  obtained 
for  the  net  present  worth  of  the  operational  costs  for  several 
time  horizons  and  interest  rates.  The  respective  present  worth 
calculations  are  performed  using  the  method  of  Gaussian  Quadrature 
for  the  numerical  integrations. 


CHAPTER  III 


SYSTEM  CONFIGURATION  AND  OVERALL  COST  MODEL 
System  Configuration 

The  system  configuration  is  depicted  in  Figure  3.1.  The  system 
consists  of  N identical  units,  each  with  P component  types.  Each  of 
the  P component  types  is  supported  with  an  inventory  of  r«j  spares, 

1 <j  <P. 

Component  Type  1 2 j ...  P 


N Units 


m nio  m.  nip 

o o o o 

o o o 
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Figure  3.1  General  System  Configuration 


The  schematic  illustration  of  Figure  3.1  could  represent  a 
lleul  of  aircraft,  cars,  or  ships,  in  which  the  component  types 
represent  basic  subdivisions  of  the  entire  vehicle,  or  unit.  However, 
the  schematic  is  a more  practical  and  realistic  depiction  of  a 
group  of  identical  radio  receivers,  each  composed  of  P "plug-in" 
modules.  Obviously,  any  particular  class  of  radios  has  a fixed 
number  of  modules  required  for  proper  operation.  This  fixed  number 
of  required  modules  is  P. 

For  a system  of  radios,  for  example,  it  is  reasonable  for  each 

of  the  P major  modules,  or  components  as  they  will  now  be  called, 

to  be  supported  by  an  inventory  of  spares.  There  are  initially 
til  ! 

spares  for  the  j component  type,  1 <j  <P.  Qualitatively,  it  would 
be  assumed  that  those  components  with  the  largest  failure  rates 
would  be  stocked  with  the  largest  number  of  spares;  however,  this 
may  not  be  necessarily  true,  since  system  availability,  or  the 
probability  that  a specified  number  of  units  is  operational  at  a 
specific  time,  is  not  the  major  criterion  for  determining  the  best 
system  configuration.  The  present  worth  of  the  system's  operational 
costs  is  the  major  criterion  of  system  worth. 

As  was  specified  in  the  introduction,  the  cost  model  for  study 
in  this  research  is  comprised  of  the  sum  of  the  present  worths  of 
the  manufacturing  or  procurement  costs,  the  repair  or  replacement 
costs,  and  the  cannibalization  costs;  the  system  revenue  or  return 
function  is  considered  as  a negative  cost  to  offset  the  other  costs. 
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Manufacturing  Costs 

The  manufacturing,  or  procurement  .costs  are  assumed  to 
be  a linear  function  of  the  number  of  components  of  each  type. 

Thus,  if  a • is  the  unit  cost  of  producing  or  procuring  the  j 

J 

component  type,  then  the  cost  incurred  at  time  zero  is: 

(3.1)  K , = N .?  a,  + J a,m.  . 

I j=l  J j~ I J J 

K-j  is  the  present  worth  of  the  manufacturing  costs  incurred  in 
order  to  begin  system  operation.  The  first  term  represents  the 
cost  incurred  to  obtain  the  working  system  of  il  independent  and 
identical  units,  while  the  second  expression  is  the  cost  of  the 
spares  for  all  of  the  P component  types. 

The  expression  given  previously  for  the  initial  manufacturing 
costs.  Equation  (3.1),  neglects  such  factors  as  quantity  discounts, 
in  which  the  component  costs  of  the  j component  type,  a.,  decreases 
for  large  quantities  procured.  It  also  neglects  any  production 
"learning  curve"  effects  in  which  the  component  production  cost 
decreases  with  increases  in  the  production  yield;  the  increased 
production  yield  is  partially  a result  of  improved  manufacturing 
techniques. 

In  applying  this  model  for  the  manufacturing  costs  it  is 
possible  to  calculate  in  such  a way  that  it  incorporates  all  the 
fixed  and  variable  costs  of  manufacturing,  as  well  as  the  trans- 
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portation,  set-up,  and  any  other  costs  incurred  upon  initially 
deploying  the  system. 

On  a present  worth  basis,  the  manufacturing  costs  will  be  valued 
at  their  initial  level  because  they  are  incurred  only  at  time 
zero.  It  is  assumed  that  no  additional  spares  are  manufactured 
or  procured  after  the  initial  system  is  set  up.  Thus,  the  present 
worth  of  the  manufacturing  costs  does  not  depend  upon  the 
nominal  annual  interest  rate,  r,  nor  the  time  horizon,  T.  Equation 
(3.1)  is  the  expression  for  the  present  worth  of  the  initial 
cost  of  manufacturing  or  procuring  the  system  depicted  in  Figure  3.1, 

Replacement  Costs 

The  replacement  costs  are  the  costs  incurred  with  replacing 
failed  components  with  components  from  the  respective  component 
spare  inventory.  Each  component  type,  ,i , will  have  associated 

with  it  a replacement  cost,  b.(t),  which,  in  general,  may  be  a 

3 

function  of  time.  The  coefficient  b^(t)  is  the  component 
replacement  cost  for  a single  replacement  of  the  j^1  component 
type.  It  is  primarily  a personnel  or  labor  intensive  cost,  since 
its  greatest  portion  is  assumed  to  be  composed  of  the  wages 
and  incidental  personnel  costs  of  those  performing  the  replacement. 

The  replacement  cost  is  defined  as  the  cost  incurred  in  locating 
the  failed  component  within  the  failed  unit,  removing  it,  and 
installing  a new  spare.  The  removed  component  is  not  repaired, 
but  discarded;  it  is  assumed  to  have  zero  salvage  value  in  the 
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s 


failed  state. 

Tlie  total  component  replacement  cost  is  a function  of 
both  the  cost  to  replace  a particular  component  type,  and  the  proba- 
bility that  the  specific  component  lias  been  repaired.  In  the  context 
of  this  analysis,  the  terms  "repair"  and  "replacement"  will  be 
used  synonymously. 

In  the  case  of  a finite  repair  time  governed  by  the  probability 

density  function  g.(t),  the  total  expected  cost  .of  repair  for 
J 

the  system  configuration  of  Figure  3.1  is: 

P N+m-j-1  ( i, repairs  of  the 

(3.2)  E(Repair  Cost(t))  = E b.(t)  E i Pr^4  jtn  component  type: 

j=l  J i=l  ’^by  time  t 

The  coefficient  b^(t)  is  the  cost  of  a single  replacement 

of  the  j ^ component  type. 

Figure  3.2  contains  the  chronological  description  of  one  failure- 
repair  cycle  occurring  before  time  t.  The  probability  of  such 
an  occurrence  for  the  /^component  type  is  calculated  as  follows: 


failure 


repair  no  failures 


L. 

0 


v 

J JL_ 

x x+dx 


y y+dy  t 


Figure  3.2  Chronological  Description  of  Exactly  1 Failure-Repair 
Completion  by  Time  t. 
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- Failure  between  "i 

x and  x+dx,  and  repair  j /-,\  (1) 

(3.3)  Pr  between  y and  y+dy,  and  = f.  ' (x)dxg  . (y-x)dyR  ■ (t-y), 

no  failures  between  ; J J J 

I y and  t.  J 

Integrating  over  all  0<x<_y,  and  0 < y < t,  one  obtains 

\'\  failure]  t v (1)  , (l), 

(3.4)  Pr  J and  > = / r f,  (x)  g.(y-x)  R-  ^ (t-y)  d xdy 

|l  repair  J 0 0 J J J 

= f^  ^(t)  * g (t)  * R.^  ^(t),  where 
j j J 

f-^(t)  is  the  probability  density  function  of  the  time  to  the 
J 

first  failure,  and  R.^(t)  is  the  effective  reliability  function 

J ^ 

associated  with  f.^(t).  The  replacement  probability  is  a 
J 

multiple  convolution,  where  the  convolution  operation  is  defined 
as  follows: 


(3.5)  A ( t ) * B(t) 


/ A(x)  B(t-x)  dx,  t>  0. 
o 


In  general,  convolution  is  commutative,  distributive,  and  associ- 
ative as  illustrated  by  the  following  equations,  which  are  presented 
for  reference  purposes  without  proof: 

(3.6)  A(t)  * C(t)  = B(t)  * A(t)  . 


(3.7)  A(t)  * ( 13 ( t ) + C(t))  - A(t)  * 0(t)  + A(t)  * C(t) 


(3.8)  A(t)  * (B(t)  * C(t))  = (A(t)  * D(t) ) * C(t)  . 


^ Appendix  I contains  the  derivation  of  f^^(t)  and  R / ^(t) 
as  used  in  this  analysis. 


- — u ■ 


i i i.MI 
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Continuing,  the  probability  of  exactly  2 failures  and  2 
repairs  for  the  component  type,  by  time  t,  is: 


;2  failures! 
(3.9)  Pr  ' and 

! 2 repairs 


= f,(1)(t)*gi(t)*f.(1)(t)*g.(t)*R.(1)(t) 
J J J J J 

= (fj(1)(t)*gj(t))*2  * Rj(1)(t). 


Figure  3.3  contains  the  description  of  all  possible  failure  and 
repair  modes  for  the  j*’1'  component  type.  Allowing  cannibalization, 
there  can  be  as  many  as  N+m--l  failure-repair  cycles  for  the 

J 

component  type,  versus  m.  replacements  without  cannibalization. 

J 

In  general,  the  probability  of  exactly  i failures,  i repair 
completions,  with  the  component  being  operational  by  time  t is: 

(1)  *i  (1) 

(3.10)  Pr  {i  failures,  i repairs}=  (f.  ( t )*g .. ( t ) ) * R . v (t), 

J J j 

(1) 

where  f.  (t)  and  g.(t)  are  the  component's  ordered  failure  and 
J J 

repair  density  functions,  respectively;  R.  (t)  is  the  "ordered" 

th  J 

reliability  function  of  the  j component  type.  The  superscript, 

■'  th 

*i , indicates  the  i multiple  convolution. 

Substituting  the  expression  for  the  probability  of  i failures 
and  i repair  completions  into  Equation  (3.2)  yields  a closed  form 
expression  for  the  expected  repair  cost  as  a function  of  time. 

(3.11)  E(Repair  Cost(t))  = y b.(t)  ^i(f  .(t)  * g.(t))  * R.(t). 

jf-1  J i=l  J J j 


In  general,  the  present  worth  of  a system's  repair/ replacement 
cost  is  obtained  by  integrating  the  product  of  the  expected  repair 


Figure  3.3  Description  of  Possible  Failure-Repair  Cycles  for  the  j Component  Type. 
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cost,  with  the  term  exp(-rt)  over  the  required  time  horizon,  M . 
The  term  exp(-rt)  is  the  continuous  annual  compounding  discount 
factor,  where  r is  the  nominal  annual  interest  rate. 

The  present  worth  of  the  system's  repair/replacement  costs 
is  designated  K^. 


(3.12)  K2  = 
T p 

/ Y b.(t) 
o j=1  J 


N+Zj'[{  i(f.(1)(t)*g-(t))*1}  * R (1)(t)]exp(-rt)  dt. 
1=1  J J J 


Cannibalization  Costs 

The  total  expected  cost  of  cannibalization  for  a system  composed 
of  N units,  depends  upon  the  probability  of  cannibalizing  a unit 
as  well  as  the  cannibalization  cost  incurred  when  a particular 
component  fails  causing  cannibalization.  The  cannibalization 
probability  is  a function  of  the  number  of  units,  N,  and  the  number 
of  spares,  mj  , for  each  component  type,  j. 

The  total  expected  cannibalization  cost  is  expressed  as  a 
sum  of  products  of  the  component  cannibalization  cost,  c j , and 
the  probability  of  that  particular  component  failing  and  causing 
cannibalization  of  unit  i,  D. .(t)„  There  can  be  as  many  as  (N-l) 

1 J 

cannibal izations  of  any  system  initially  composed  of  fl  units. 

N-l  V 

(3.13)  E(Canni bnl ization  Cost  (t))  = Y Y Cj  D ..  (t) 

i=l  j=l  J 1J 
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Since  each  unit  is  assumed  identical,  and  the  periodic  replace- 
ment of  failed  components  does  not  significantly  alt'"''  the  config- 
uration of  any  unit,  then  the  cannibalization  probability  is 
the  same  for  any  unit  within  the  system.  It  only  varies  with 

the  component  type,  j,  1 <j<P . Thus  D. .(t)  = D.(t).  The  cannibal- 

J J 

ization  probability,  D^t)  will  be  expressed  using  the  density 
functions  of  the  ordered  times  to  failure  of  the  ^component  type. 
The  revised  cannibalization  cost  equation  is: 

P 

(3.14)  E(Cannibal ization  Cost(t))  = (N-l)  £ c.D.(t) 

j=l  J J 

This  revised  equation  results  from  the  basic  premise  that  as  many 
as  (N-l)  units  may  be  cannibalized,  each  with  equal  1 i kl i hood. 

Stochastic  Description  of  Cannibalization 

The  probability  expression,  D.(t),  will  be  expressed  using  the 

vJ 

marginal  density  functions  of  the  ordered  failure  times  of  the 
component  type,  1 <j_<P . It  can  be  shown,  however,  that  for  a very 
large  time  horizon,  the  probability  of  cannibalizing  a unit  because 
of  failure  of  a particular  component  type  is  not  a function  of 
time,  but  depends  upon  the  number  of  units,  N,  the  number  of  spares, 
m • , and  the  respective  hazard  rate  for  each  component  type  j,  l<j<P. 

J 

For  a system  initially  composed  of  N units,  with  P components 
per  unit,  arid  m. spares  for  the  j^‘  component  type,  l<j£P,  a unit 
is  cannibalized  if  (N+m--!)  failures  occur  in  component  type  j 
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before  (N+rn,-  1)  failures  occur  in  any  other  component  type,  k, 


with  k/j,  l£k<P.  Thus,  D.(t)  is  the  probability  that  the  (N+m.-l) 

th  ^ ^ 

failure  time  of  the  j component  type  is  less  than  the  smallest 


th 


of  the  (N+m^-1)  failure  times  of  the  remaining  component  types; 


then  it  is,  by  definition,  less  than  all  of  the  (N+rn^-1)  failure- 


times  of  the  remaining  (P-1)  components. 


The  following  analysis  is  based  upon  independent  identically 


th 


distributed  times  to  failure  for  the  j component  type.  Define 


.th 


t (M+nij  1)  as  the  tl-me  t0  the  (N+m  i)  failure  of  the  j 

^ th 

component  type,  and  D.(t)  as  the  probability  that  the  j compo- 

J 


nent  type  caused  the  cannibalization. 


(N+m.-l ) 

(3.15)  D . (t)  = Pr  { t.  J < 
J <3 


(N+m.-l ) 

Min(tk  k ,.,.,)} 


k=l ,. . ,P 
k/j 


If  the  (N+m.-l  )^  failure  of  the  j^  component  type  is  less 

J 


than  the  minimum  of  the  corresponding  failure  times  of  the 
remaining  components,  then  it  is  less  than  all  of  the  corresponding 
failure  times  of  the  remaining  components. 


(3.16)  D (t)  = Pr  { t J < All ( t k ),.,.,)} 

j j k 

k=l ,. . ,P 

k/j 


Under  the  assumption  of  independent  component  failures,  the 


, .th  th 

probability  that  the  (N+m.-l)  failure  time  of  the  j component 

J 


type  is  less  than  all  of  the  corresponding  failure  times  of  the 
remaining  component  types  is  the  product  of  the  probabilities 


iv 


- 
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that  the  (N+m.-l)**1  failure  time  of  the  j*'*1  component  type  is  less 

J 

than  each  of  the  (N+m.-l)  failure  times  of  the  k component  types,  for 

k = 1 , 2,  ... , P,  k i j. 


(3.17)  D.(t)  = n Pr  {t.(N+mj'^  < t^N+mk_1h  . 

J k=1  J k 

kj^j 

The  probability  expression  given  in  Equation  (3.17)  is  calculat- 


ed using  the  marginal  density  function  of  the  time  to  the  (N+m.-l) 

J 

failure  of  the  j*'*1  component  type.  The  required  marginal  density 


th 


function  of  the  ordered  failure  times  is  obtained  by  integrating  the 
joint  density  function  of  the  respective  component  ordered  failure 
times;  the  joint  density  function  is 

h(t-| ( N+,t>i -1 ) ^ ^(N+n^-l) > tp^,Hmp"^).  Because  of  the  inde- 

pendence of  the  respective  ordered  failure  times,  the  joint  density 
function  of  the  ordered  times  to  failure  is  the  product  of  the  re- 
spective marginal  density  functions,  h^*4|71l"^(t),  h2^+m2~^(t),..., 
etc. 


Since  the  marginal  density  functions  of  the  ordered  times  to 
failure  of  each  of  the  component  types  are  proper  density  functions, 
the  definite  integrals  of  the  marginal  density  functions, 
h^  (^+ml-l ) (t) , h2^+IT12"^  (t) , etc.,  from  zero  to  infinity,  all  re- 
duce to  unity.  Consequently,  the  following  double  integration 
results: 


r 


(3.18)  Pr  { t.^N+inj'^  < t (^"V1)  } _ 

J k 

ft  fy  h (N+mj-1 ) (x)  h (N+mk-l) 
o o J k 


(y)  dx  dy. 


Using  previous  reasoning,  the  following  expression  yields  the  time 
dependent  probability  of  a cannibalization  due  to  a failure  of  the 

x L. 

jun  component  type. 


(3.19)  D.(t)  = n Z1 /y  h.(N+mj“1)(x)  hk(N+V1)(y)  dxdy  . 

k/j 

If  the  mean  times  to  repair  a component  type  are  significantly 
less  than  the  respective  mean  time  between  failures,  it  is  possible 
to  make  use  of  an  instantaneous  replenishment  model,  in  which  the 
marginal  density  functions,  h.^(t),  for  j = 1,  2,  P,  and 

J 

i = 1,  2,  ....N+m.-l ..,  are  equivalent  to  the  probability  density 
functions  of  the  i**1  ordered  failure  of  the  j^1  component  type,  i.e. 

fj(i)(t). 


Case  of  Finite  Replacement  Time 


In  the  case  of  a finite  replacement  time  governed  by  the  density 

function  g.(t),  the  joint  density  function  of  the  ordered  times  to 
J th 

the  (N+m.-l)  failure  must  include  the  time  to  complete  the  replace- 
ment of  each  failed  component  type.  The  marginal  density  function 


of  the  time  to  the  i^1  failure  of  the  component  type,  h.^(t), 

J 

is  the  density  function  of  the  sum  of  the  times  of  the  i ordered 
failures,  and  the  (i-1)  repairs  of  the  jth  component  type. 

Under  the  condition  of  a non-zero  replacement  time,  or  finite 
replacement  rate,  with  replacement  time  governed  by  the  density 
function  g-(t),  it  is  not  necessarily  true  that  h.^(t)  = f.^(t). 

J ' J J 

The  density  function  of  the  sum  of  the  i independent  failure  times 
and  (i-1)  independent  repair  times  is  the  convolution  of  the  i-fuld 
convolution  of  the  ordered  failure  density  with  the  (i-l)-fold 
convolution  of  the  repair  density. 


(3.20) 


(i) 


(t)  = 


J 


g.(t) 

J 


;(i-D 


l<j<P, 


l<i<M+m.. 

J 


For  example,  the  time  to  the  third  failure  of  the  first  compo- 
nent type  will  be  the  sum  of  the  first  three  failure  times  and  the 
first  two  repair  times  of  the  first  component  type.  In  this  case, 
the  marginal  density  function  of  the  time  to  the  third  failure  of 
the  first  component  type  is  expressed  as  a multiple  convolution  of 
the  ordered  failure  and  repair  density  functions. 


(3.21)  ' h/3)(t)  = f1U/(t)*g1(t)*f/U(t)*g1(t)*f1U/(t) 

(1)  *3  *? 

• Y (t)  * 9i(t)  • 


(i) 


o), 


o), 


In  terms  of  the  marginal  densities  of  the  terminal  failure 
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times,  h . (|sj+mj ~1 ) ( t)  and  h.  ^ (t) , the  cannibal  ization  proba- 

J * 

bility,  D-(t),  is  expressed  as  follows: 

J 

(3.22)  D . (t)  = n f1  fy  h (N+mj_1)(x)  h^^V^ty)  dxdy  . 

J k=J  O u J K 

k£i 

Substituting  the  correct  form  of  the  convolved  failure  and  repair 
density  functions  into  the  marginal  density  functions  yields  the 
correct  expression  for  the  cannibalization  probability  function: 


(3.23) 


Dj(t) 


P 

n 

k=l 

k^j 


o o 


m *(N+mi-l ) *(N+mr2) 

(f,U,(x)  J * g.(x)  J ) X 

J 3 


(fk(1,(y)‘(N^-1) 


, *(N+m,,-2l 

(y)  k ) dx  dy. 


Present  Worth  of  Cannibalization  Cost 

Define  K3  as  the  present  worth  of  the  cannibalization  cost. 

Assume  c.  is  the  cost  of  cannibalizing  a unit  due  to  failure  of 
<3 

the  jth  component  type,  which  occurs  with  probability  Dj(t).  Also 
assume  that  the  nominal  annual  interest  rate  is  r,  and  that  interest 
is  compounded  continuously. 

Based  upon  Equation  (3.14),  the  present  worth  of  the  expected 
cannibalization  costs,  over  a time  horizon  T,  is  expressed  as 
follows: 
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T P 

(3.24)  « = / (N-1 ) Z c,D-(t)  exp(-rt)  dt„ 

3 0 j = l J J 

is  explicitly  a function  of  the  system  size,  as  measured  by 

the  number  of  units,  N,  and  the  number  of  components  per  unit, 

P.  The  cost  is  implicitly  a function  of  the  spare  inventory  size, 

m.  , since  the  term  D.(t)  is  a function  of  both  N,  and  m.. 

J J J 

The  coefficient  c.  is  determined  empirically  from  previous 
J 

cannibalization  activity  on  a given  system.  The  time  horizon  is 
estimated  by  the  system  user,  and  will  be  based  upon  previous 
data  on  similar  systems.  It  could  also  be  determined  by  the  system 
designer  before  the  specification  requirements  for  the  system  are 
written. 

System  Revenue  or  Return  Functions 

The  basis  for  calculating  a system's  return  or  worth  to  the 
user  is  the  present  worth  of  the  instantaneous  revenue  function. 

The  system  revenue,  or  return  function  V.j(t),  is  a monetary  measure 
of  the  value  of  i units  to  the  system  user  when  they  are  operational 
at  time  t.  If  it  is  not  measured  in  "dollars  per  unit  time",  for 
example  a production  process  whose  output  is  measured  in  "output 
per  unit  time",  then  it  should  be  converted  to  a monetary  measure 
before  using  this  model. 

Brown f7’  justified  a monotonically  decreasing  revenue  function 
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with  increasing  age,  on  the  basis  that  increasing  obsolescence  and 
decreasing  efficiency  will  usually  cause  a monotonic  decrease  in 
the  revenue  function.  Also  with  the  attrition  of  units,  the  revenue 
function  of  the  system  usually  decreases,  and  continued  depletion 
of  the  system  could  result  in  a system  having  too  few  units  to  operate 
successfully.  This  is  exemplified  by  a Taxi  'Service  having  too  many 
vehicles  undergoing  repair;  the  result  would  be  long  waiting  times 
for  current  customers  who  may  utilmately  switch  to  another  Taxi 
system,  resulting  in  decreased  revenue  for  the  original  Taxi  service. 

Even  though  the  system's  instantaneous  revenue  function,  V..(t), 
will  usually  increase  monotonically  for  any  increase  in  the  number 
of  operating  units,  and  the  revenue  function  could  increase  pro- 
portionately for  initial  increases  in  the  number  of  units,  additional 
units  will  eventually  be  of  less  value. 

The  law  of  diminishing  marginal  returns  states  that  as  equal 
increments  of  one  input  are  added  to  a production  system,  with  the 
inputs  to  the  other  productive  services  (remaining  units)  being 
held  constant,  the  resulting  incremental  production  (system  revenue) 
v/i  1 1 decrease.  However,  the  law  of  diminishing  marginal  returns 
says  nothing  about  the  effects  of  increasing  the  number  of  units 
whenever  technological  processes,  which  could  alter  the  revenue 
producing  capability  of  subsequent  units,  are  changing.  The  units 
are  assumed  to  be  operating  in  a stable  technological  environment 
throughout  the  time  period  of  interest,  so  the  law  of  diminishing 


marginal  returns  is  applicable. 

A consequence  of  the  lav;  of  diminishing  marginal  returns  is 
to  partially  limit,  or  provide  an  upper  bound  to,  the  number  of 
units  initially  deployed,  N;  the  upper  bound  is  the  number  of  units 
for  which  the  marginal  operating  cost  of  the  system  just  equals 
the  marginal  revenue  derived  from  system  operation.  This  micro- 
economic  result  would  fix  the  optimal  number  of  units  initially 
deployed  to  a point  where  the  marginal  revenue  obtained  from  op- 
erating the  Nth  Unit  was  equal  to  the  marginal  cost  of  operating 
that  unit.  The  marginal  cost  would  be  the  incremental  manufactur- 
ing costs,  plus  the  expected  repair  costs  for  the  last  unit,  and 
the  expected  cannibalization  costs  for  the  last  unit  deployed. 

However,  this  microeconomic  analysis  does  not  account  for  the 
stochastic  system  behavior.  Specifically,  it  assumes  that  the 
availability  of  each  of  the  units  is  unity;  the  system  would  be 
assumed  to  operate  deterministically,  with  N units  available  through 
out  its  life. 

Since  the  unit  availability  function  P (t),  is  the  probability 

7 

of  exactly  i units  being  operational  at  time  t,  the  revenue  obtain- 
ed from  having  i units  operational  is  simply  the  product  of  the 

unit  revenue  function,  V (t),  and  the  unit  availability  function, 

i 

P.(t)  , i.e.  Vi(t)  Pitt). 

The  unit  revenue  function,  V\(t),  is  said  to  be  a monotonically 
increasing  function  of  the  number  of  units,  i,  if  for  any  two  values 


i'l  and  i'2,  with  i^  < the  relation  V 1 ( t ) < V^(t)  holds.  V..(t) 
would  be  strictly  monotonically  increasing  if  only  the  strict 
inequality  holds  true.  Conversely,  V^(t)  is  monotonically  decreas- 
ing if  and  only  if  -V- (t)  is  monotonically  increasing. 

It  is  assumed  that  since  the  time  dependency  of  the  system 
revenue  is  accounted  for  with  the  unit  availability  function,  then 
the  revenue  function,  V.j(t),  can  be  represented  as  a time  indepen- 
dent monotonically  increasing  function  of  the  number  of  units.  Thus, 

(3.25)  Vi(t)  = Vi. 

However,  the  specific  form  of  the  revenue  function  will  depend 
upon  the  particular  system  and  its  usage.  The  system  user  is 
expected  to  obtain  the  specific  revenue  function  which  is  appro- 
priate. This  can  be  done  either  analytically,  for  some  systems,  or 
empirically,  for  others. 

The  present  worth  of  the  expected  revenue  obtained  over  a 
time  duration  T,  is  designated  K^. 

T N 

(3.26)  K4  = / E V-  P.:(t)  exp(-rt)  dt. 

0 i-0  ( 1 

Ttie  summation  extends  over  the  range  of  units,  from  0 to  M, 
because  even  though  at  least  1 unit  is  usually  required  for 
a positive  return,  there  may  be  some  cases  where  a negative 
return  results  when  0 units  are  operational.  The  model  is 
applicable  in  such  cases  by  choosing  negative. 


The  calculation  of  P^t),  the  probability  of  exactly  i units 
operating  at  time  t,  depends  upon  the  specific  spare  configuration 
considered  and  will  be  left  for  the  example  problem  of  Chapter  V. 

Example  Calculations 

Consider  the  schematic  illustration  of  Figure  3.4,  in  which 
N=2,  P=3,  rn-|=2,  m2=l , and  m3=l . 


Lb 

□ 

1 

LJ- 

-0-0 

0 

o 

c 

c 

Figure  3.4  Specific  System  Configuration 

The  expected  manufacturing  cost  is  equivalent  to  its  present 
worth  since  the  units  and  their  spare  supplies  are  manufactured  at 
time  zero.  Thus  the  present  worth  of  the  expected  manufacturing 
costs  is  given  by: 

3 3 

(3.27)  K,  = 2 l a.  + l a-m. 

1 j=l  J j=l  J J 

= (2+mj)a-|  + (2+m2)a2  + (2+013)33 


= 4a i + 3a0  + 3a 


For  a fixed  system  configuration,  the  present  worth  of  the  total 
manufacturing  or  procurement  costs  is  a linear  function  of  the 
respective  component  procurement  costs. 


The  expected  repair  cost  for  the  system  shown  in  Figure  3.4 
(p.  40)  is  calculated  below  using  Equation  (3.11). 

3 1+m 


(3.28)  E( Repair  Cost(t))  = I b,  i(f  (1  *(t)*g  .(t))*1  * R.(1)(t) 

j=l  Ji=l  J J J 

= b^f  (1)*9l*R  (1)  + 2(f  (1)*9])*2*R i(1) 

/n  /n 


+ 3(f1(1)*g1f3*R1(1)} 
+b2{f2(l)*g2*R2(1)  + 2(f2(1)*g2)*2*R2(1)) 

+b3{f3(1)*g3*R3(1)  +2(f3{1)*g3)"2*R3(1)} 


For  notational  convenience,  the  time  dependence  of  the  failure, 
repair,  and  reliability  density  functions  was  not  given  in 
Equation  (3.28). 

The  present  worth  of  the  expected  repair  cost  for  the  system 
of  Figure  3.4,  with  interest  compounded  continuously  at  a nominal 
annual  rate  r,  over  a time  horizon  T,  is  given  below. 

(3.29)  Ko  = / l bi1rj[i(f.(1)(t)*g.(t))*1  * R {1)(t)]oxp(-rt)  dt 
'•  0 j-1  Ji=l  J J J 


Assuming  an  exponential  density  function  of  the  time  to  failure 
and  time  to  repair,  it  is  possible  to  obtain  the  closed  form  expres- 
sion for  the  present  worth  of  the  expected  repair/replacement  costs 


for  the  system  of  Figure  3.4(p.  40). 

The  specific  density  functions  of  the  time  to  failure,  time 
to  repair,  as  well  as  the  reliability  function  are  given  below. 


(3.30) 

f.(,) 

J 

(3.31) 

g.(t) 

(3.32) 

R(1) 

J 

f.  (t)  = 3,  exp(-Bit),  where 
J 0 J 

3.  = NX-,  for  j=l , 2,  and  3. 
J J 


Using  the  basic  formulas  of  Table  3.1,  the  following  convo- 
lution results  are  obtained. 


(3.33)  f.0)(t)  * g.(t)  = — (exp(-p.t)  - exp(-B.t)} 

j j 3j  ‘ J J 


} t exp(-Xjt) , ii  Bj  - Uj ■ 


(11  *2 
(f.U)(t)  * g (t)) 

J J 


P • p . 

( 8-  Zt~}  ( t(exp(-u-t)  + exp(-B- 1) ) 
Pj  Mj  J J 

-ifeSPL-Hj t ) _-_expX- fiji ]) } 

(Bj  - Vj) 


(t'/B)  exp(-y-t),  if  8-  = P-. 

J J J 


r 


1 


n ^ *3 

(3.35)  (f.  (t)*g,(t)) 

J *> 


= (.A5)'i-)3((t2/2)(exp(-pit)  - cxp(-S.t)) 

V “j  J ! 

- (3exp(-v.t)  + 3exp(-B.t)) 

(Bj-liJ  3 3 

6 y (exp(-u,t)  - exp(-B.t))} 

T Q J J 

= (pj6/3! ) (t5/20)  exp(-Pjt) , if  3j=Pj. 


The  respective  reliability  convolutions  are  given  below. 


(3.36,  fW(,|w  - (&>< 

J j j Uj  j 

- t exp(-6,t)} 

J 

= yi2(t2/2)  exp ( -vi  t) , if  3.=Vj • 

J J J 

m *2  (i)  4 

(3.37)  (f.  *g  .(t))  *R  (t)  = 

J J J 

/e.u.f  ^pxnf-K  t)  - 3exp(-p.t)  + 2t_ex[^Ji.jtV  t | cx|)(-p,j,t) 


.Wl 


S - "j 


Bd  - fj- 


+ (t/2)  exp(-Sjt)  J, 


p,2(t2  /2)  exp(-pjt),  if  p «B.. 
! J J 


r 


45 


+ (1/2)  ( t(exp(-u  t)  - 3t2  exp(-B.t))  - t3  expf-0 .t) 1 

J('0j  - Mj)  12  ? 


The  present  worth  of  the  expected  repair/repl acement  costs 
is  obtained  by  integrating  the  product  of  the  expected  repair 
costs  and  exp(-rt),  where  r is  the  continuous  annual  interest 
rate.  Equation  (3.29)  contains  such  a calculation  for  the  system 
of  Figure  3.4  (p.  40) . 

Using  Equation  (3.28),  together  with  the  results  derived 
in  Equations  (3.36-3.38),  the  expected  repair  cost  for  the  system  of 
Figure  3.4  is  tabulated  in  Tables  3.2  and  3.3.  Table  3.3  contains 
the  expected  repair  cost  equation  when  the  system  hazard  rate, 

6.  (=  MX.),  equals  the  repair  rate,  p . , for  j=l , 2,  and  3. 

J J I j 

Assuming  b^-l^b^l,  the  normalized  present  worth  of  the 
repair/rcplaccmcnt  costs  is  calculated  below.  A comparison  will  be 


n 
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made  among  the  present  worth  values  for  particular  hazard  and  repair 
rates.  In  the  case  of  exponential  failure  and  repair  densities,  the 
analysis  will  be  tabulated  as  a function  of  the  system  mean  time 
between  failure,  and  the  mean  time  to  repair,  MTBF  and  MTTR,  respec- 
tively. The  system  mean  time  between  failure,  for  the  component 

type,  is  1/B-,  which  equals  1/NX.;.  Also,  assuming  eouivalent  hazard 
J y 

and  repair  rates  for  each  component  type,  the  present  worth  of  the 
expected  repair/replacement  costs  is  calculated'  using  the  equations 
given  in  Table  3.3  (p.  49). 

From  Equation  (3.29),  the  closed  form  expression  for  ^ is: 

T 3 

(3.39)  y =/  b Z i(f1(1)(t)*g1(t))*i  * R,(1)(t)  exp(-rt)  dt 
o 1 i=l  1 1 1 

T 2 (11  *i  (1) 

+'  b2  2 i(f'  ;(t)*gJt))  * R_  (t)  exp(-rt)  dt 

o i=l  ^ ^ 

T 2 (l)  *i  (l) 

+/  bo  Z 1(f  (t)*go(t))  * R (t)  exp(-rt)  dt. 

o 3i=l  3 3 3 

Substituting  the  formulas  of  Table  3.3  into  Equation  (3.29), 
and  normalizing  the  replacement  cost  coefficients,  b.=l,  for  j=l , 2, 
and  3,  yields  the  following  integral  equation  for  the  present  worth 
of  the  repair/replacement  costs. 
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1 ? o 

(3.40}  = / y1  (r  /2)  expC-^  + r)t)  dt 


T 2 2 

+ / y„  (t  /2)  exp(-(y  + r)t)  dt 
o ^ 

T 2 2 

+ / y3  (t  12)  exp (-(u3+  r)t)  dt 

o ° 

T 4 4 

+ 2 (/  y1  (t  /6'4)  exp(-(y1  + r)t)  dt 

o 

T 4 4 

+ 2 (/  y (t  /6-4  ) exp(-(y  + r)t)  dt 

o 2 ^ 

T 4 4 

+ 2 (/  u3  (t  /6’4  ) exp(-(u3  + r)t)  dt 
o 

T 6 6 

+ 3 / y-j  (t  / 20 - 6 * 3 ! ) exp(-(y1  + r)t)  dt 
o 


Tables  3.4  and  3.5  contain  listings  of  the  normalized  present 
worth  of  the  repair/ replacement  costs,  Kg,  for  the  system  of  Figure 
3.4,  as  a function  of  time,  and  the  continuous  interest  rate.  Annual 
time  periods  require  the  dimensions  of  the  interest  rate  to  be 
"percent  per  year".  By  "normalized  present  worth",  is  meant  the 
present  worth  with  cost  coefficients,  b-j,  bg,  and  b^,  of  unity. 

Table  3.4  contains  the  values  of  K?  with  the  MTBF  and  MTTR  both  equal 

i c 

to  1.00,  while  Table  3.5  contains  the  values  when  both  the  MTBF, 
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and  MTTR  equal  10.00  . They  were  calculated  using  Equation  (3.40) 


Examining  the  tables  indicates  that  the  normalized  present 
worth  of  the  repair/replacement  costs  is  relatively  insensitive  to 
variations  in  the  interest  rate  between  10%  and  20%  per  period.  Also 
noted  is  the  fact  that,  for  specified  time  horizons  between  1 and 
10  periods,  the  normalized  present  worth  under  the  larger  values 
of  MTBF  and  MTTR  is  significantly  less  than  that  obtained  with  the 
smaller  values  of  MTBF  and  MTTR.  The  reason  for  this  is  that  with  the 
larger  values  of  MTBF  and  MTTR,  there  is  a smaller  probability  of 
failing  and  completing  the  repair  or  replacement  procedure  within 


the  specified  time  horizon.  Consequently,  the  expected  cost  of  the 
failure  and  repair  completion  is  small,  as  is  its  present  worth. 

Example  Calculation  of  Cannibalization  Probability 

Consider  again  the  system  shown  in  Figure  3.4  (p.  40),  with 
2 units,  3 components  per  unit,  and  2 spares  for  the  first  com- 
ponent type,  1 spare  for  the  second,  and  1 spare  for  the  third. 

According  to  the  original  definition,  D^t)  is  the  probability 
of  cannibalizing  a unit  because  of  a failure  of  the  first  component 
type.  In  the  example  of  Figure  3.4,  a unit  is  cannibalized  due  to 
a failure  of  type  1 only  if  the  third  failure  of  component  type  1 
occurs  before  the  smaller  of  the  second  failure  times  of  components 
2 and  3.  Specifically, 

2 

Appendix  II  contains  tables  for  the  present  worth  calculations 
of  l<2,  for  a system  MTBF  of  1.00,  and  an  MTTR  between  0.10  and  0.33. 
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(3.41)  D i ( t ) = Pr  { t/3^  < Min(t^^2^  , t^2)  ) } . 

(3.42)  D (t)  = Pr  { t2(2)  < MinU^3)  , t^2)  ) } . 

(3.43)  D3(t)  = Pr  { t^2)  < Min(t^3^  , t ^ ) } . 


( n\ 

Here,  t (3)  is  the  third  failure  time  of  component  type  1,  t2 

1 (2) 
is  the  second  failure  time  of  component  type  2,  and  t3  is  the 

second  failure  time  of  component  type  3.  The  probability  that  a par- 
ticular order  statistic  is  less  than  the  smaller  of  two  other  order 
statistics  is  the  probability  that  it  is  smaller  than  both  of  the 

other  order  statistics.  Thus  the  expressions  for  D.(t),  for  j=l, 

J 

2,  and  3,  are  as  follows: 


(3.44)  D1  (t)  « Pr  (t^3^  < t^2\  and  t^3^^  } 

(3.45)  D2(t)  = Pr  {t  2{Z)  <t1(3)  , and  t2(2)<  t3(2)} 

(3.46)  D3(t)  = Pr  {t3^2^  < t^3^  , and  t3^<  t^2)} 

Since  the  components  are  assumed  to  fail  independently  according 

to  the  specific  failure  density  function,  f.(t),  then  any  order  sta- 

J 

tistics  of  different  component  types  must  be  independent.  Consequent- 
ly. the  joint  events  expressed  in  equations  (3.44-3.46)  are 


r 


Independent. 


(3.47)  D-|  (t)  = Pr  {t/3)  < t2(2))  X Pr  [t/3)  < t3(2)) 


(3.48)  , D2(t)  = Pr  (t2(2)  < X Pr  (t^^  < t^^} 

(2)  (3)  (2)  (2) 

(3.49)  D3(t)  = Pr  {t3  < t1  } X Pr { t3  < t2  } 


Using  the  assumption  of  independence  among  each  of  the  respec- 

(3)  (2)  (2) 

tive  ordered  times  to  failure,  t^  , t2  , and  t3  , the  joint 
density  function  of  the  ordered  failure  times  is  the  product  of  the 
marginal  densities.  The  cannibalization  probability  expressions  given 
in  Equations  (3.47-3.49)  are  expressed  as  the  product  of  the  integrals 
of  the  joint  density  functions.' 


(3.50)  D]  (t)  = / h2(2,(y)  dx  dy  X 


f f h^3^(x)  h ^(y)  dx  dy. 
o o * 3 


(3.51)  D„(t)  = / /\(2)(x)  h.(3)(y)  dx  dy  X 

c.  0 0 2 1 


t y (2)  (?) 

I f h2  (x)  hg  (y)  dx  dy. 


X 


(3.52)  D,(t)  = f r IV  ;(x)  h/  '(y)  dx  dy 

j o o J 1 

t y (2)  (2) 

/ / h.U'(x)  h„Uj(y)  dx  dy. 

0 0 J L 

The  above  results  are  consistent  with  the  general  expression 

for  D.(t)  given  in  Equation  (3.19),  even  though  they  were  derived 
J 

independently. 

For  the  case  of  an  infinite  replenishment  rate,  the  marginal 
density  functions  of  the  ordered  time  to  failure  for  each  component 
type  are  equivalent  to  the  respective  ordered  failure  time  density 
functions,  f-j^(t),  f 2^^ (t ) , and 

(3.53)  h^3)(t)  = f^3,(t). 

(2)  (2) 

(3.54)  h2  (t)  = f2  (t). 

(2)  (2) 

(3.55)  h3  (t)  = f3V  '(t). 

The  equations  for  the  density  function  of  the  ordered  time  to 
failure  are  given  below. 


(3.5G) 

fl(3,(t>  - 

( ^ ^ 

F^t)2 

^(t) 

(1-F-j  (t) ) . 

(3.57) 

! 2 A 

t . . „ ! 

F?(t) 

f9(t) 

(1-F2(t))°. 

r 


(3.58)  f3^(t)  = 


1,1,0  / 


F3(t)  f3(t)  (1-F3(t))  . 


For  the  exponential  density  function  of  time  to  failure,  with 

the  hazard  rate  for  each  component  type,  j,  given  by  the  marginal 

J 

density  functions  of  ordered  times  to  failure  are: 


(3)  2 

(3.59)  h^  ;(t)  = 12  X1  (l-exp(-X  t))  exp(-2X  t) 


(3.60)  h V '(t)  = 2 \9  (l-exp(-.\„t) ) exp(-lA0t) 


(3.61)  h3  (t)  = 2 A3  (l-exp(-X  t))  exp(-lX  t) 


Substituting  the  expressions  given  in  Equations  (3.59-3.61) 

into  the  equations  for  the  component  cannibalization  probability, 

I 

(3.50-3.52),  and  performing  the  respective  integrations,  yields  the 
time  dependent  component  cannibalization  probability,  D-(t),  for 

J 

j=l , 2,  and  3.  These  expressions  for  the  time  dependent  cannibaliza- 
tion probability  are  plotted  in  Figures  3.5  through  3.10,  as  a 
function  of  time.  Note  that  they  do  not  integrate  to  one  because 
they  are  not  plotted  as  probability  density  functions;  D.(t)  repre- 

J 

' sents  the  time  dependent  probability  that  the  component  type 
caused  a cannibalization.  Each  set  of  graphs  contains  the  respective 
cannibalization  probability  function  for  selected  values  of  the 
MTCF  for  the  system  configuration  of  Figure  3.4  (p.  40). 


Component  Cannibalization  Probability,  D (t),  versus 
MTBF(1 )=MTBF(2)=MTBF(3)=1 .0.  J 


Component  Cannibalization  Probability,  D.(t),  versus 
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By  varying  the  magnitude  of  the  respective  failure  rates,  (re- 
ciprocal MTBF's)  the  time  required  to  reach  a steady  state  probabil- 
ity level  varies.  This  is  seen  by  comparing  Figures  3.5  and  3.10. 

In  Figure  3.5  (p.  59),  each  component  type  has  an  MTBF  of  1.0  and 
the  steady  state  probability  level  is  reached  within  approximately 
3 time  periods.  In  Figure  3.10  (p.  64),  however,  the  respective 
component  MTBF  values  are  all  2.0  and  the  steady  state  probability 
level  is  reached  after  6 time  periods. 

The  steady  state  cannibalization  probability  is  obtained  by 
evaluating  the  cannibalization  probability  equations  with  the  upper 
limit  of  the  outer  integrals  set  at  infinity. 

Table  3.6  contains  selected  steady  state  cannibalization 

probability  values  as  a function  of  the  respective  component  hazard 

rates,  for  the  system  configuration  of  Figure  3.4,  having  N=2,  P=3, 

m-|=2,  ni2=l , and  1113  = ! . It  also  contains  the  probability  of  no 

cannibalization,  1 -D-j -D2-D3 . The  tabulated  values  indicate  that  it 

is  primarily  the  relative  magnitude  of  the  hazard  rate,  rather  than 

the  absolute  value,  which  determines  the  respective  cannibalization 
3 

probability.  For  example,  when  the  respective  hazard  rates  are 
all  increased  or  decreased  by  the  same  relative  magnitude,  the 
corresponding  cannibalization  probability  does  not  change.  This 
is  seen  when  all  of  the  hazard  rates  are  equal  to  1 , or  100. 


Appendix  III  contains  the  steady  state  probability  derivation. 
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Discussion  and  Conclusions 

This  chapter  has  outlined  the  stochastic  description  of  cannibal- 
ization as  a maintainability  technique  and  contains  a specific  oper- 
ational cost  model.  Only  those  costs  which  were  directly  related 
to  the  basic  system  configuration  were  considered;  specifically 
included  were  the  manufacturing  or  procurement  costs,  the  repair 
or  replacement  costs,  and  the  cannibalization  costs.  A revenue  or 
return  function  was  included  in  the  model  as  a negative  cost. 

An  example  calculation  was  performed  for  a specific  system 
configuration  and  each  of  the  costs  were  analyzed  separately  on  a 
present  worth  basis.  Plots  of  the  cannibalization  probability  as 
a function  of  time  were  given  and  its  limiting  steady  state  behavior 
was  found  to  be  only  a function  of  the  relative  magnitude  of  the 
respective  component  hazard  rates. 

Also  included  was  a discussion  of  the  important  factors  considered 
in  the  system  revenue/return  function.  No  specific  analytic 
expression  was  given  for  the  revenue/return  function  because  the 
return  to  the  system  user  is  a function  of  the  specific  application 
required,  as  we 1 1 as  the  usage  environment  encountered. 

Chapter  IV  contains  a discussion  of  specific  system  constraints 
and  the  example  calculations  required  for  a constraint  on  spare 
adequacy. 
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CHAPTER  IV 
SYSTEM  CONSTRAINTS 


Introduction 

Since  sparing  allocation  is  an  integer  process,  i.e.  there  can 
only  be  an  integer  number  of  spares  allocated  to  any  component  type, 
the  classical  methods  of  constrained  optimization,  such  as  the 
Jacobian  method  and  the  Lagrangian  method  [43, 50],  are  not  applicable 
since  they  assume  a continuous  domain  of  spare  component  types. 

Under  the  formulation  specified  in  Chapter  I,  in  order  for 
a specific  spare  allocation  to  be  optimal,  it  must  result  in  a 
minimum  total  operational  cost  and,  in  some  cases,  satisfy  one  or 
more  constraints.  If  the  constraints  are  applicable,  they  are  usually 
imposed  by  the  user  during  the  system's  design  stage.  The 
constraints  are  usually  in  the  form  of  performance  requirements,  or 
limitations  on  available  resources. 

Performance  constraints  are  generally,  but  not  necessarily, 
formulated  using  inequalities  of  the  form  "greater  than  or  equal 
to  (>0",  while  resource  constraints  are  usually  formulated  using 
"less  than  or  equal  to  (<)"  inequalities. 

Typical  performance  constraints  are  minimum  requirements  on  a 
system's  reliability  over  a specific  time  horizon,  or  a requirement 
on  its  interval  availability  over  the  same  time  horizon.  In  this 


69 


context,  the  reliability  is  defined  as  the  probability  that  the  system 
operates  over  some  pre-selected  time  interval,  while  the  interval 
availability,  as  defined  by  Rau  [34],  is  the  expected  fraction  of 
an  interval  of  specified  length  that  the  system  is  in  an  operable 
state.  Another  performance  constraint,  which  will  be  discussed  in 
detail,  is  a constraint  on  the  probability  of  spare  adequacy  over 
a specified  time  horizon.  This  constraint  requires  that  the  probabil- 
ity of  experiencing  a critical  spare  shortage  before  some  specified 
time  be  strictly  less  than  a specified  value,  e.g.,  l-aQ.  By 
"critical  spare  shortage"  is  meant  an  inventory  depletion  which 
precludes  further  system  operation.  Obviously,  this  performance 
constraint  can  be  expressed  conversely  as  a requirement  that  the 
probability  of  experiencing  no  critical  spare  shortages  be  at  least 
as  large  as  aQ,  the  complement  of  the  previously  specified  constraint 
value. 

Some  obvious  resource  constraints  are  limitations  on  the  total 
weight  or  volume  of  the  spare  components;  another  is  a limitation 
on  the  total  dollar  investment  in  spares.  The  resource  constraints 
appropriate  to  this  sparing  allocation  problem  can  be  quantitatively 
described,  in  general,  as  follows: 

Assume  the  maximum  allowable  consumption  of  the  k*^  resource, 

e.g.  weight,  volume,  or  available  funds,  is  Q^,  then  the  feasible 

solution  space  of  spare  component  allocations,  (m^ , m^,  ...,  nip), 

L d k j (mj ) < Qk,  k - 1 ,2,..S}. 

J 1 


belongs  to  the  set  L,  where  L = {nij : 
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The  function  ^ kj (mj ) yields  the  quantity  of  the  ktfl  resource  consumed 

X L. 

by  the  allocation  of  mj  spares  for  the  jLr  component  type,  and  the 
index  k varies  from  1 to  S for  constraints  on  S total  resources. 

Under  the  required  resource  constraint,  the  problem  is  to  find 
the  allocation  of  spares,  (m^ , m^,  ...,  mp),  which  minimizes  the 
appropriate  cost  function,  and  belongs  to  the  set  L. 

Until  now,  nothing  has  been  said  about  the  effects  of  a finite 
time  horizon  on  the  optimal  spare  allocation.  Obviously,  the  length 
of  time  the  system  is  expected  to  be  operating  is  significant  in 
determining  the  optimal  spare  allocation.  Certain  types  of  constraints 
are  not  meaningful  for  certain  applications;  for  example,  in  the 
analysis  of  a tactical  weapons  system,  the  time  horizon  may  be 
very  short,  consequently,  the  use  of  a steady  state  availability 
constraint  is  obviously  inappropriate;  in  fact  it  is  meaningless 
unless  an  infinite  number  of  replacements  are  permitted  for  a 
particular  class  of  components,  requiring  an  allocation  of  an 
infinite  number  of  spares  for  that  component  type. 

It  became  necessary  to  develop  a realistic  constraint  on  the 
spares  which  was  responsive  to  changes  in  the  system's  expected 
usage  life,  T,  as  well  as  to  changes  in  the  number  of  fielded  units, 

N.  An  appropriate  performance  constraint  applicable  over  a finite 
time  horizon  is  a constraint  on  the  probability  of  experiencing  a 
critical  spare  shortage  for  any  of  the  component  types. 
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Constraint  on  Spare  Shortages 

It  is  necessary  to  obtain  the  minimum  number  of  spares  required 
in  order  to  meet  a specific  constraint  on  the  probability  that  there 
will  be  no  critical  spare  shortage  for  any  of  the  P component  types. 
One  such  constraint  is  the  requirement  that  the  probability  that  no 
critical  spare  shortages  occur  for  any  of  the  P component  types  be 
at  least  as  large  as  some  user  specified  value,  say  a0-  A critical 
spare  shortage  was  defined  previously  as  an  inventory  c .pletion 
which  precludes  further  system  operation.  Specifically,  system 
operation  terminates  when  more  than  (N+m--l)  failures  of  the  jth 

J 

component  type  occur,  for  j = 1,  2,  ...,  P.  Under  cannibalization, 
as  many  as  (N-l)  components  may  be  added  to  the  spare  inventory  of 
any  one  component  type. 

Mathematically,  this  constraint  is  defined  as  follows: 

)no  critical  spare  shortage  | 

is  experienced  for  any  of  the  P component  > _>  aQ  . 

types  during  the  interval  [0,TJ  j 

Define  PL  as  the  maximum  number  of  possible  replacements  of 

J 

the  component  type.  If  cannibalization  is  permitted,  Mj  = N+nij-1  , 
otherwise  Mj  = mj.  Assuming  component  failures  occur  in  accordance 
with  the  Poisson  distribution,  with  mean  time  between  failures,  1/Aj, 
and  assuming  a finite  time  horizon,  T,  the  probability  of  experiencing 
more  than  failures  of  the  component  type  is  no  greater  than 
the  probability  of  experiencing  more  than  nij  failures,  with  N>  1. 
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This  implies  that  it  is  easier  to  protect  against  critical  spare 
shortages  when  cannibalization  is  permitted  than  when  it  is  not, 
since  the  system  can  tolerate  as  many  as  (N+m.-l)  failures  rather 

J 

iL 

than  just  mj  failures  for  the  j component  type. 

It  may  seem  superfluous  to  impose  a constraint  on  spare 
adequacy,  when  the  cannibalization  process  itself  is  used  as  a tool 
for  preventing  critical  spare  shortages.  However,  the  definition  of 
a critical  spare  shortage  was  a spare  inventory  depletion  which 
prevented  further  system  operation.  Under  cannibalization,  the  point 
at  which  system  operation  terminates  is  extended  because  of  the 
permissible  use  of  cannibalized  components  to  repair  impending 
component  failures. 

For  any  given  component  type,  j,  the  probability  that  no  more 
than  Mj  spares  are  required  during  the  interval  [o,t]  is  just  the 
probability  that  the  (Mj+l)^  failure  time  exceeds  T.  If  h j ^ f ^ ^ ( t ) 

1.  L 

is  the  marginal  density  function  of  the  time  to  the  failure  of 
the  j**1  component  type,  then  the  probability  that  Mj  spares  are 
adequate  is  just  the  probability  that  there  will  be  no  more  than 
Mj  failures  in  [o,t] . Alternatively,  this  is  just  the  probability 
that  the  (Mj+l)th  failure  time  exceeds  T.  Thus  the  constraint  is 
written  as  follows: 

+h  00  (Mt+1  ) 

(4.2)  Pr  {(M-j+l)th  failure  time  > T)  = / h-  J (t)  dt. 
j j j 

With  P statistically  independent  components  per  unit,  it  is 
necessary  that  the  (Mj+l)tfl  failure  time  exceed  T for  all  component 


types.  Thus  the  constraint  becomes: 

(4.3)  II  / h.  J U(t)  dt  > a . 

j=l  T J “ 

If  there  are  N independent  and  identical  units,  the  effective 
hazard  rate  for  the  j**1  component  type  is  NX.,  where  Xi  is  the 
hazard  rate  of  a single  component  of  the  j™  type. 

Assuming  Poisson  failures,  the  times  between  failures  are 
exponentially  distributed.  With  exponential  failure  times  and 

instantaneous  replacement,  the  probability  density  function  of  the 
sum  of  (Mj+1)  independent  and  identically  distributed  random 
variables,  h-^'j+^(t),  is  the  gamma  density  function. 

J 

(4.4)  hj%fl)(t)  = NXj(NXjt)Mj  exp(-NXjt)  / Mj ! 

The  instantaneous  replacement  assumption  is  justified  if  the  mean 
time  to  replace  a particular  component  type  is  significantly  less 
than  its  mean  time  between  failures,  under  the  instantaneous 
replacement  assumption,  it  is  not  necessary  to  include  the  times 
required  to  perforin  the  replacements  in  the  marginal  density  function 
of  the  time  to  the  (M,+l)th  failure  of  the  jth  type. 

J 

Integrating  Equation  (4.4)  by  parts  yields: 

CO  [ij_. 

(4.5)  / (NAJ  (NAit)Mj  exp(-NX^t)  dt  = IJ  (NXJ;k  exp(-NAJ). 

T J J k=0  J k 1 

M j ! K U K- 


Substituting  the  result  of  Equation  (4.5)  into  Equation  (4.3)  yields 
the  complete  constraint. 


(4.6)  n Ej  (NXiT)k  exp(-NX -T)  > a. 

j=l  k=0  J rH- 


where  aQ  is  the  minimum  probability  of  spare  adequacy,  as  specified 

by  the  system  user,  and  Mj=rnj  when  cannibalization  is  not  permitted, 

or  M.=N+m--l,  when  it  is. 

J ^ 

The  level  of  protection  against  critical  spare  shortages  for 
the  component  type  is  the  cumulative  Poisson  sum  given  in 
Equation  (4.5).  The  constraint,  as  written  in-Equation  (4.6)  is  the 
product  of  independent  cumulative  Poisson  densities,  each  with 


\ 


mean,  Nx-T. 

sj 

If  one  does  not  impose  the  requirement  that  any  one  component 
type  has  a higher  priority  for  protection  against  critical  spare 
shortages  than  other  component  types,  and  if  there  are  no  restric- 
tions on  the  maximum  number  of  spares  al 1 owed  for  any  component 
type,  then  one  can  allocate  spares  for  each  of  the  P component  types 
in  accordance  with  Equation  (4.6)  as  follows: 

1.  Starting  with  the  first  component  type,  solve 
for  M-j  by  setting  M3 , . . . , Mp  = °°.  This  forces 
the  Poisson  summation  over  components  2 
through  P to  unity,  and  one  need  only 


N 


* 


find  the  smallest  value  of  M-j  which  satisfies: 

Mi  k 

E (li’X.T)  exp(-KXiT)  > a . 
k=0  1 k!  ' 0 
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This  can  be  calculated  directly,  or  by  using  a 
table  of  the  Cumulative  Poisson  Distribution 
Function. 

2.  Next,  determine  the  smallest  value  of  Mo  which 
satisfies: 

m2  k 

(4.8)  Z (NAoT)  exp(-KA9T)  > a . 

k=0  6 ~k!  0 

3.  Continue  this  procedure  for  the  remaining  com- 
ponents. The  minimum  number  of  spares  provisioned 

t h 

for  the  Pu  component  type  will  satisfy: 

Mp  k 

(4.9)  Z (tUpT)  exp(-NAnT)  > an. 

k=0  k!  r ~ 0 

This  procedure  is  used  to  calculate  a lower  bound  on  the 
required  number  of  spares  for  each  component  type.  The  allocation 
achieved  will  not  necessarily  satisfy  the  constraint  of  Equation 
(4.6),  because  the  maximum  of  each  Poisson  summation  is  1.0  when 
M.  equals  infinity;  thus  the  product  over  all  P component  types 

vJ 

will  be  at  most  1.0.  (For  almost  all  applications,  it  will  be 
significantly  less  than  1.0.)  However,  if  an  exhaustive  search  is  to 
be  performed  over  the  total  costs  associated  with  each  allocation, 

then  this  procedure  will  yield  the  starting  allocations  over  which 
the  search  should  be  conducted.  Then,  for  each  allocation,  the 
product  of  the  respecti ve^Poisson  sums  can  he  calculated,  in  accor- 
dance with  Equation  (4.6),  to  eliminate  those  allocations  which  do 


r 


not  satisfy  the  constraint  given  in  Equation  (4.6). 

The  outlined  procedure  can  be  performed  in  any  order,  i.e. 
starting  with  any  component  type.  However,  it  is  only  applicable 
if  no  component  type  has  a higher  priority  for  spare  adequacy  than 
other  component  types. 

If  one  particular  component  type,  say  type  "q",  requires  a 

spare  sufficiency  level  of  a for  example,  then  M should  be 

q M 

calculated  as  the  smallest  integer  which  satisfies: 


m . 

(4.10)  r1  (MATT  exp(-NAT)  > ac 


The  remaining  minimum  spare  allocations  are  calculated  in  accordance 
with  Equations  (4.7)  through  (4.9). 

If  all  of  the  components  do  not  have  distinct  MTBF’s,  for 
example  if  two  or  more  of  the  components  have  equivalent  MTBF's, 
then  they  will  necessarily  have  identical  Poisson  mean  values,  and 
will  be  allocated  the  same  number  of  spares  initially. 

There  may  be  some  situations  in  which  the  spare  adequacy 
constraint  given  in  Equation  (4.1)  is  further  restricted  with  the 
requirement  that  each  of  the  component  types  must  have  the  same 
level  of  protection  against  critical  spare  shortages . With  such  a 
restriction,  the  constraint  given  in  Equation  (4.6)  becomes: 

M-  k 1/p 

(4.11)  T. 3 (NX  .T)  exp(-NA.T)  > (a..)  , for  j = 1,  ...,  P, 

i^-n  J "IT  At  — u 


where  M.=in.,  when  cannibalization  is  not  permitted,  and  H . - ( M f m • - 1 ) , 
J J J »■' 

when  it  is.  Thus,  the  spare  adequacy  probability  is  partitioned 


V. 


equally  among  the  P component  types. 


The  equal  partitioning  is  justified  as  long  as  no  component 

type  has  a higher  priority  for  spare  sufficiency  than  any  other 

component  type.  If  there  were  an  additional  requirement  that  the 

t h 

q component  type  have  a spare  adequacy  probability  in  excess  of 

a , then  the  single  constraint  of  Equation  (4.11)  becomes  the 
H 

multiple  constraint: 


(4.12) 


Zq  (NA  T)k  exp (-NX  T) 

k=0  q k!  q 


M • 

£**  (NA  T)k  exp  (-NX -T)  > for  j = 1 , . . . ,P 

k=0  J k!  ~J“  - d q j i q. 


Tables  4.1  through  4.5  contain  the  minimum  spare  allocations 
in  order  to  meet  a minimum  probability  of  spare  adequacy, a , in 
excess  of  0.70,  0.75,  0.80,  0.85,  and  0.90,  for  various  values  of 
N,  MTBF,  and  T.  Each  table  lists  the  minimum  spare  allocations 
necessary  under  both  conditions:  when  cannibalization  is  permitted, 
and  when  it  is  not.  The  values  tabulated  are  the  minimum  number 
of  spares  required,  itk  , for  each  component  type,  j,  j = 1,  2,  and  3, 
calculated  in  accordance  with  Equations  (4.7)  through  (4.9). 

The  spare  allocations  are  the  same  in  Tables  4.4  and  4.5,  because 
the  respective  values  of  NX -T  are  the  same  in  both  tables. 

The  tabular  solutions  listed  are  relevant  in  this  application 
because  they  are  dependent  upon  both  the  time  horizon,  and 
the  failure  rates  of  each  of  the  respective  component  types. 


TA3LE  4.1  Minimum  Spare  Allocations  for  a 2 Unit  System  with  a Time  Horizon  of  5 Periods. 
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Constraint  on  Total  Weight  or  Volume  of  Spares 


If  it  is  desired  to  analyze  the  operational  cost  model  subject 
to  a constraint  on  the  total  weight  or  volume  of  the  spare  inventory, 
the  formulation  would  be  similar  to  the  classical  "knapsack  problem" 
of  integer  programming  [47],  However,  in  the  "knapsack  problem"  the 
value  to  the  user  is  a linear  function  of  the  essentiality,  or  worth, 
of  each  item.  This  is  not  necessarily  the  case  in  the  problem  pre- 
sented in  this  dissertation;  the  present  worth  of  the  revenue  func- 
tion is  a nonlinear  function  of  the  number  of  units. 

A constraint  on  the  total  weight  of  the  spare  inventory  would 
be  applicable  in  a situation  such  as  modular  sparing  of  processing 
units  for  a computer  on-board  an  aircraft.  Here,  both  the  total 
weight  and  total  volume  may  have  restrictions  applied  to  them. 

If  both  constraints  are  in  effect,  the  required  spare  allocation 
would  have  to  satisfy: 


(4.13)  £ w,m,  < W0,  and 

j=l  J J 


f , V;  i v 


Hero,  w.  is  the  unit  weight  of  the  j11  component  type,  UQ  is  the 
maximum  weight  allowed  for  all  P component  types,  v.  is  the  volume 
of  the  j 1 component  type,  and  VQ  is  the  maximum  allowable  volume 
for  all  P component  types. 
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j 

Constraint  on  Total  Investment  in  Spares 

j 

In  some  economic  situations,  there  may  be  a budget  limitation 
which  restricts  the  available  funds  for  purchasing  the  spare  inven- 
tory components. 

This  constraint  may  be  total,  in  that  no  more  than  L dollars 

P 

may  be  spent  for  the  entire  spare  inventory  of  z m-  components, 

j=1  3 

or  it  may  be  partial,  in  that  a maximum  amount  may  be  spent  only  on 
"critical"  spare  components,  with  no  funds  allocated  for  other 
"less  critical"  component  types.  This  restriction  limits  the  spare 
inventory  allocation  to  only  those  components  which  have  available 
funding. 

Under  the  first  situation,  in  which  the  total  inventory  dollar 
value  is  limited  to  Lq  dollars,  the  appropriate  constraint  is: 

P 

(4.14)  z 

j=l  3 3 

Here,  a.  is  the  cost  of  the  jth  component  type,  excluding  any 
J 

discounts  for  large  quantities  purchased. 

If  there  is  the  restriction  that  only  certain  spare  types  may 
be  purchased  with  the  limited  total  budget  of  LQ  dollars,  then  the 
spare  allocation  scheme  is  more  complex.  Specifically,  if  only  spare 
component  type  "q"  may  be  purchased  because  of  its  importance,  or 
its  criticality  for  mission  success,  then  the  appropriate  constraint 
i s : 
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(4-15)  V’qlLo- 

This  specifics  a limit  on  the  number  of  spares  of  component  type  q, 

(4.16)  niq  < L0/aq, 

where  aq  is  the  cost  of  a single  component  of  type  q. 

If  only  two  component  types  are  permitted,  say  q and  r,  the 
appropriate  constraint  is: 

(4.17)  aqmq  + arinr  < L0, 

with  m , and  mr  restricted  to  integers. 

Constraint  on  Interval  Availability 


Since  component  reliability  does  not  allow  consideration  for 
component  repair  or  replacement  actions,  and  system  reliability 
is  not  an  adequate  indication  of  system  dependability  when  repairs 
at  the  system  level  are  permissible  it  is  necessary  to  develop  a 
meaningful  constraint  on  a measure  of  system  dependability  which 
allows  for  the  possibility  of  component  repair.  Interval 
Availability  is  such  a measure. 

Interval  Availability  was  defined  previously  as  the  expected 


fractional  amount  of  an  interval  of  specified  length  that  the  system 
is  in  an  operable  state. 


Given  a specified  time  horizon,  T,  over  which  the  system  is 
expected  to  be  operating,  a constraint  can  be  imposed  on  the  cost 
model  requiring  that  the  system's  (all  N units)  interval  availability 
be  in  excess  of  aQ,  for  example.  Thus  it  is  required  that  the 
expected  fraction  of  the  time  horizon,  T,  in  which  all  N units 
are  operable,  be  in  excess  of  aQ. 

However,  in  order  to  calculate  the  interval  availability,  one 
needs  an  accurate  measure  of  the  time  dependent  system  availability, 
A(t).  This  may  be  difficult  to  formulate  analytically  for  many 
system  configurations.  The  interval  availability,  Aj(T),  is  the 
average  value  of  the  system  availability  over  a finite  time  hor- 
izon, T. 

T 

(4.18)  At(T)  = (1/T)  / A(t)  dt. 

0 

Thus,  an  appropriate  constraint  on  interval  availability  is: 

(4.19)  Aj(T)>a0. 


Maintainability  Constraint 


According  to  Ran  [34],  system  maintainability  is  a character- 
istic of  design  and  installation  which  is  expressed  as  the  proba- 
bility that  an  item  will  be  restored  to  specified  conditions  1 i' 
a given  period  of  time  when  maintenance  actions  are  perfor 
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Under  some  circumstances,  for  example,  the  repair  of  emer- 
gency generating  equipment  within  a hospital,  there  may  be  a 
constraint  imposed  on  the  maximum  allowable  downtime.  One 
formulation  of  this  maintainability  constraint  is: 

. 

(4.20)  Pr  (downtime  j<td}  _>  a0, 

where  aQ  is  the  minimum  required  probability  of  experiencing  a 
downtime  no  greater  than  t^. 

System  downtime  is  usually  defined  as  the  total  time  from 
system  failure  until  return  to  service  and  includes  the  time  to 
perform  the  replacement;  or  repair,  as  well  as  any  administrative 
time;  the  administrative  portion  of  the  downtime  is  that  required 
to  procure  spares  and  process  any  manuals  or  test  equipment  nec- 
essary for  the  repair  or  replacement  action. 

The  major  element  of  downtime  which  may  be  decreased  by  having 
an  adequate  supply  of  spares  (in  accordance  with  the  criterion 
set  forth  in  Equation  (4.6))  is  the  portion  of  the  administrative 
downtime  which  includes  the  lead  time  required  when  an  additional 
spare  procurement  request  is  being  processed. 

However,  if  the  administrative  time  is  relatively  small  in 
comparison  with  the  replacement  time,  then  the  system  downtime 
can  be  approximated  by  the  component  replacement  time.  In  the 
context  of  this  research,  the  administrative  portion  of  downtime 
will  bo  assumed  small  relative  to  the  replacement  time,  thus  the 
total  downtime  and  total  replacement  time  are  essentially 


equivalent.  Consequently,  the  chosen  spare  allocation  does  not 
significantly  affect  the  expected  downtime  for  a system  under  this 
formulation.  Much  more  significant  are  the  statistical  parameters 
of  the  density  function  of  the  time  to  perform  the  replacements. 

Summary  and  Discussion 

This  chapter  outlined  some  specific  constraints  which  should 
be  considered  during  a system's  design  formulation.  Several  con- 
straints were  discussed,  and  a constraint  on  the  probability  of 
spare  adequacy  was  analyzed  in  detail. 

The  constraint  on  spare  adequacy  was  realistic  in  this  applica- 
tion since  the  process  of  cannibalization  is  directly  dependent  upon, 
and  in  turn  affects,  the  spare  inventory  level  of  each  of  the  com- 
ponent types.  The  use  of  cannibalization,  as  illustrated  by  the 
tables,  permitted  a system  configuration  with  fewer  total  spares  for 
each  component  type. 

Other  constraints  mentioned  were  a constraint  on  the  total 
weight  or  volume  of  a specific  spare  allocation,  a resource  constraint 
on  the  total  investment  in  spares,  a requirement  on  a system's 
interval  availability,  and  a maintainability  constraint  on  system 
downtime.  Each  of  these  constraints  may  be  applicable  to  a specific 
system  under  varying  circumstances.  The  system  user  must  decide 
initially,  which  if  any,  of  the  constraints,  or  combinations  of 
constraints  are  applicable  to  his  situation. 
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CHAPTER  V 
SOLUTION 


Model  Requirements 

In  order  to  effectively  utilize  this  model,  one  needs  specific 
cost  and  performance  information  on  the  system.  It  is  necessary  to 
know  the  component  cost  parameters;  these  may  be  obtained  from 
prior  purchase  invoices  or,  on  new  systems,  from  vendor  price 
estimates  for  a particular  purchase  quantity  at  a given  time.  One 
must  also  specify  the  number  of  units,  N,  which  are  to  be  fielded  as 
a complete  system,  as  well  as  estimates  of  the  unit  revenue  func- 
tion, , i.e.,  the  incremental  revenue  obtained  when  i units 
are  in  operation.  If  a system  does  not  generate  income,  then  a 
monetary  measure  of  its  operational  worth  to  the  user  should  be 
estimated.  This  estimate  would  then  be  used  in  the  corresponding 
calculation  of  the  present  worth  of  the  system's  return  or  revenue 
function,  K^. 

As  far  as  the  performance  characteristics  of  the  system  are 
concerned,  it  is  necessary  to  obtain  estimates  of  the  underlying 
failure  and  replacement  density  functions,  f^(t),  and  gj(t),  re- 
spectively, as  well  as  their  statistical  parameters;  the  density 
function  of  the  time  to  the  first  failure,  f.^  ^(t)  is  obtained 

J 


from  f.(t)  in  accordance  with  the  procedure  outlined  in  Appendix  I. 

These  functions  could  be  obtained  easily  enough  on  currently  fielded 

I 

systems,  by  setting  up  an  efficient  data  base  collection  system; 

1 

however,  for  a new  system  such  information  is  minimal,  and  one  would 
have  to  use  sound  engineering  judgment  to  estimate  the  system's 
performance  characteristics.  This  is  the  major  drawback  in  the  use 
of  this  model  on  new  systems. 

Another  requirement  is  that  the  net  cost  of  alternative  spare 
allocations  must  be  compared  over  the  same  time  horizon,  T;  this 
is  because  the  cost  model  is  a present  v/orth  technique. 

As  discussed  in  Chapter  IV,  there  may  be  several  constraints, 
or  limitations,  on  the  spare  inventory  allocation  which  greatly 
restrict  the  number  of  alternative  allocations  which  must  be  ex- 
amined. It  is  generally  advisable  to  impose  the  constraint  re- 
quirements first,  in  order  to  determine  the  feasible  set  of  spare 
allocations,  before  performing  any  specific  cost  calculations.  This 
will  greatly  reduce  the  total  number  of  calculations  and  comparisons 
required.  However,  this  method  is  only  applicable  if  the  active 
constraints  limit  the  feasible  spare  allocations  to  a finite  and 
comparatively  small  set.  Otherwise,  one  would  have  to  obtain  an 
upper  bound  on  the  maximum  number  of  spares  for  each  component 
type,  if  possible,  and  perform  a search  over  the  feasible  region. 

If  it  were  not  possible  to  obtain  this  upper  bound,  one  would 

I 

have  to  iterate  the  complete  cost  model  across  successive  spare 
allocations  and  determine  the  trend  in  the  cost  calculations. 
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Nevertheless,  this  technique  would  only  lead  to  a locally  optimal 
solution,  with  no  guarantee  that  the  achieved  solution  was  a global 
optimum. 

A specific  example  problem  will  now  be  discussed  in  detail  to 
illustrate  the  application  of  this  particular  cost  model.  There  are 
two  constraints  which  limit  the  feasible  spare  allocations  to  four 
possibilities.  A constraint  on  the  minimum  spare  adequacy  proba- 
bility determines  a lower  bound  on  the  number  of  spares  for  each 
component  type,  while  an  investment  constraint  on  the  total  dollar 
value  of  the  spares  provides  an  upper  limit  to  the  permissible 
number  of  spare  allocations. 

Example  Problem 

Consider  the  case  of  a small  computer  peripheral  equipment 
leasing  firm.  Available  Computer  Enterprises  (ACE).  The  firm 
operates  a long  term  leasing  service  in  which  a support  subsystem 
of  computer  peripherals  is  leased  for  an  annual  charge  of  $8000.00; 
the  basic  subsystem  consists  of  a High  Speed  Line  Printer,  a Drum 
Plotter,  and  an  Optical  Card  Reader. 

ACE,  being  in  the  embrionic  stage  of  its  development  in  the 
computer  leasing  industry,  has  only  two  systems  currently  leased. 
Besides  supplying  the  basic  peripheral  subsystems,  ACE  is  respon- 
sible for  all  maintenance  on  the  fielded  systems,  i.e.  repairing 
or  replacing  failed  peripheral  components. 


ACE  has  also  included  in  the  leasing  contract  a warranty 
provision  stipulating  that  the  customer's  risk  of  experiencing  a 
shutdown  due  to  a critical  spare  shortage  will  be  no  greater  than 
10%,  during  the  first  year's  service,  i.e.,  both  customers  are 
assured  that  there  is  at  least  a 90%  probability  that  there  will 
be  no  critical  spare  shortage.  Also  stipulated  in  the  contract  is 
the  provision  that  the  customers  are  not  obligated  to  pay  for  any 
time  that  the  system  is  down  for  maintenance  .or  repair.  This 
requirement  causes  the  revenue  or  return  function,  K^,  to  be 
stochastic  in  nature,  in  that  no  revenue  is  generated  when  both 
systems  are  down. 

In  order  to  meet  the  constraint  on  spare  shortages,  ACE  must 
stock  enough  backup  spares;  however,  the  firm  is  restricted  in  the 
number  of  spares  it  can  stock,  in  that  it  only  has  a budget  of 
$80,000.00  for  backup  spares. 

Cannibalization  has  been  suggested  as  a possible  technique 
for  decreasing  the  number  of  backup  spares  necessary  to  satisfy  the 
spare  sufficiency  constraint,  and  meet  the  limitation  on  the 
maximum  permissible  investment  in  spares. 

The  solution  of  this  allocation  problem  will  consist  of  the 

<H>  ^ ^ ^ 

ordered  set  (m^  , m2  , m-j  ),  where  m-j  is  the  optimal  number  of 

spares  for  the  Line  Printer,  is  the  optimal  number  of  spare 

★ 

Drum  Plotters,  and  m-j  is  the  optimal  number  of  Optical  Card  Read- 
ers to  be  stocked.  The  values  of  the  continuous  annual  interest 
rate  to  be  used  in  present  worth  cost  calculations  are:  10%,  15%, 


and  20%  per  year. 


f 


Table  5.1  contains  the  respective  procurement,  repair/replace- 


ment,  and  cannibalization  cost  coefficients  for  this  system. 

I 


TABLE  5.1  Cost  Coefficients  ($) 


Procurement  Repair/Replacement  Cannibalization 

Cost  Cost  Cost 


Line  Printer  12000.00  330.00 


1850.00 


Drum  Plotter  9050.00  250.00 


2300.00 


Optical  Card  5450.00  190.00  . 4000.00 

Reader 


Using  data  collected  on  other  peripheral  devices  (based  on 
the  same  technology),  ACE  officials  were  able  to  obtain  realistic 
estimates  of  the  failure  and  repair  rates  for  the  respective 
devices.  Table  5.2  contains  the  failure  and  replacement  rate  data 
for  each  device,  given  on  both  an  hourly  and  yearly  basis,  assum- 
ing an  operating  time  of  8 hours  per  day,  6 days  per  week,  52  weeks 
per  year;  also  included  are  the  yearly  MTBF  and  MTTR  values,  which 
were  calculated  under  the  assumption  of  exponential  failure  and 
replacement  time. 


f 


* 

* -X 

* * * 

* * * 

* * * 

* * * 

* * * 

r 

* 

* 

* 

* 

* 

* 

L 

* 

* 

* 

* 

* 

* 

* 

* 

* 

1 

♦ 

* 

o 

cr» 

cr> 

* 

ft 

* 

o 

r— 

i — 

* 

* 

^* 

o 

O 

o 

* 

* 

cc 

• * 

o 

CVJ 

cvj 

* 

>’ 

* 

\— 

s-  * 

1— 

» — 

r— 

* 

■K 

>>  * 

• 

* 

• 

* 

■ 

• 

* 

s: 

^* 

o 

o 

o 

* 

* 

* 

* 

p ' 

LQ 

* 

* 

* 

O 

* 

* 

* 

* 

* 

* 

* 

* 

* 

IS* 

c u 

t/> 

•«— 

s- 

CL 

s~ 

C D 


* * 


CQ 


* * 

* * * 

* * 

CO 

o 

o 

CVJ 

o 

o 

r— 

CVJ 

CO 

o 

o 

CVJ 

CO 

CO 

r— 

o 

r-- 

O 

CO 

s- 

CD 

4- > 
CD 
E 

CO 

5- 


* 

* 

* 

* 

* 

* 

* 


************ 


C 


cO 

* 

r—  * 

o 

o 

LO 

Q_ 

* 

1 * 

CO 

CVJ 

* 

CD 

• * 

«o- 

CVJ 

r— 

cn 

* 

4-> 

s*.  * 

CVJ 

CD 

o 

c 

* 

CO 

>>* 

• 

• 

• 

•r— 

* 

Cd 

^ * 

r— 

o 

o 

4-> 

* 

* 

CO 

* 

CD 

* 

U 

* 

S- 

* 

«=J* 

*3* 

CD 

CD 

* 

3 

—H* 

1 

1 

1 

CL 

* 

T— 

r—  * 

o 

o 

o 

o 

* 

•1“ 

1 * 

r— 

f— 

r— 

* 

CO 

• * 

X 

X 

X 

4-> 

* 

U_ 

* 

o 

LO 

o 

C 

* 

sz  * 

• 

• 

• 

CD 

* 

— "X 

LO 

CVJ 

LO 

C 

o 

CL 

E 

O 

O 


CVJ 

LO 

LlJ 

_J 

CQ 


* 

* 

* 

-X 


************ 


* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 


4-> 

* 

* 

* 

3 

* 

* * 

* * * 

* i.  * 

* * * 

* * * 

CL 

* 

* 

* 

* 

.y 

* 

o 

* 

^* 

* 

o 

* 

CD 

r—  * 

* 

* 

-*-> 

1 * 

o 

CVJ 

CVJ 

* 

CD 

* 

cO 

• * 

o 

CO 

CO 

* 

r— 

* 

Cd 

S~  * 

• 

• 

• 

* 

JD 

* 

>>* 

o 

CO 

CO 

* 

CO 

* 

+-> 

— ' -X 

r— 

* 

* 

c 

* 

* 

•r— 

* 

CD 

* 

* 

1 CO 

* 

E 

* 

* 

> 

* 

CD 

- — * 

o 

CO 

CO 

* 

* 

O 

r—  * 

o 

CO 

CO 

* 

* 

cO 

1 * 

CO 

CO 

* 

* 

i — 

• * 

o 

o 

o 

* 

o 

* 

Q- 

S-  * 

o 

o 

o 

* 

<*- 

* 

CD 

-C  * 

• 

• 

• 

* 

* 

Cd 

— *x 

o 

o 

o 

* 

tn 

* 

* 

* 

* 

* 

t- 

s- 

“O 

* 

* 

* 

CD 

CD 

S- 

* 

* 

* 

4-> 

4-> 

CO 

* 

* 

* 

c 

4-> 

O 

* 

* 

* 

•I— 

o 

* 

* 

* 

s- 

r— 

r— 

* 

* 

* 

CL 

Cd 

CO 

* 

* 

* 

U CD 

* 

* 

•X 

CD 

E 

•r-  “O 

* 

* 

* 

c 

3 

•M  cO 

* 

* 

* 

•r— 

i- 

D-  CD 

* 

* 

* 

_] 

Q 

O cd 

* 

* 

* 

* 

* * 

* * * * 

* * * 

* * * 

* * * * 

* * 

< 


95 


Using  the  respective  failure  rates,  the  minimum  sparing  alloca- 
tions required  to  satisfy  the  spare  adequacy  constraint  of  90%  for 
the  first  year's  operation  may  be  calculated;  since  the  failure 
rates  were  much  smaller  than  the  replacement  rates,  it  was  per- 
missible to  use  the  Poisson  model  described  in  Chapter  IV  to  cal- 
culate the  required  sparing  allocations.  Specifically,  it  was 
necessary  to  determine  m.  for  j = 1,  2,  and  3,  such  that: 

J 

3 Mi  k 

(5.11  n ZJ  (2X.1)  exp(-2X-il)  > 0.90,  where  M.,  = 2+m--l. 

0=1  k=0  J k!  " J J ' 

The  values  given  in  Table  5.3  are  the  minimal  spare  allocations 

which  must  be  checked  against  the  spare  adequacy  constraint  given 

in  Equation  (5.1);  allocations  are  given  for  adequacy  values  of  70%, 

75%,  80%,  85%,  and  90%.  For  a 90%  requirement,  the  minimum  spare 

allocation,  assuming  cannibalize tion  is  permitted,  is: 

(5.2)  (m-j , m2,  m3)  = (4,  2,  0). 

However,  when  the  respective  Poisson  Sums  are  substituted  into 
Equation  (5.1),  the  product  is  0.8978,  which  does  not  meet  the 
90%  constraint.  The  spare  adequacy  constraint  may  be  satisfied  by 
stocking  at  least  4 Line  Printers,  2 Drum  Plotters,  and  1 Optical 
Card  Reader;  this  yields  a spare  sufficiency  of  0.9209,  at  a 
cost  of  $71,550,  which  is  $8,450  under  the  $80,000  budget  constraint, 
quantitatively  described  as  follows: 


(5.3)  Z a,m.  < $80,000.00. 
j=l  J J 
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There  are  4 feasible  allocations  which  meet  the  spare  adequacy  con- 
straint, yet  do  not  violate  the  budget  restriction:  (4,2,1),  (4,2,2), 

(4,3,0),  and  (5,2,0). 

The  cost  model  will  be  analyzed  completely  for  each  allocation, 
over  time  horizons  of  5,  10,  and  15  years. 

The  present  worth  of  the  procurement  costs  is  not  dependent  upon 
the  interest  rate  or  time  horizon,  because  it  is  incurred  at  time  0. 


(5.4) 


K,  - 2 Z a:  + £ a^m^ 
1 j=l  J j=l  J J 


"$124,550.00  for  (4,2,1) 

$130,000.00  for  (4,2,2) 

$128,150.00  for  (4,3,0) 

. $131  ,100.00  for  (5,2,0) 

The  present  worth  of  the  repair/replacement  costs  is  calculated 

in  accordance  with  Equation  (3.12)  (p.  29),  and  requires  calculating 

(1)  *i  (1) 

the  multiple  convolution,  (f.  (t)*g.(t))  * R.  (t),  in  order  to 

J j J 

obtain  the  probability  of  exactly  i failures  and  i repair  comple- 
tions, for  i = 1,  2,  N+m.-l. 

J 

Under  an  exponential  failure  density,  it  is  possible  to  con- 


vert the  above  multiple  convolution  into  the  same  form  as  the  mar- 

t h t h 

ginal  density  functions  of  the  time  to  the  i failure  of  the  j 

component  type,  h.^(t).  From  Equation  (3.20)  (p.  34), 

(i)  (1)  J *1  *i-l  (1)  (1) 

h,  (t)  = f (t)  *g,(t)  . Since  R,  (t)  = f (tj/p,,  then 

J J J,  J J vJ 


(1), 


O), 


(5.5)  (f j vt)*9j (t) ) * Rj'"(t)  = f. 


0)(t) 

(i+1) 


M+l 


* g.(t)  /b. 

J u 


= h . ( t )/3 . » 

J J 

for  j = 1,  2,  P,  and  i = 1,  2,  N+m.-l. 
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Because  of  the  complexity  of  the  expressions  for  the  multiple 
4 

convolutions  , it  was  necc-;sdry  to  perform  the  present  worth 

5 

integrations  numerically,  using  the  method  of  Gaussian  Quadrature  . 

The  calculation  of  K^,  the  present  v/orth  of  the  expected 
cannibalization  cost  is  similar  to  the  calculation  of  K2  except 
that,  in  addition,  it  is  necessary  to  obtain  products  of  double 
integrals,  in  order  to  calculate  the  time  dependent  cannibalization 
probability;  these  multiple  integrals  were  calculated  numerically, 
using  the  method  of  Gaussian  Quadrature.  The  general  form  of  the 
cannibalization  probability  was  given  in  Equation  (3.22)  (p.  35). 

After  the  time  dependent  cannibalization  probability  is 
calculated,  one  must  multiply  by  the  respective  cannibalization 
cost  coefficients,  ci5  and  the  continuous  compounding  discount 
factor,  exp(-rt),  and  then  integrate  over  the  required  time  horizon, 
T,  again  using  Gaussian  Quadrature. 

For  the  4 possible  configurations,  the  present  worth  of  the 
expected  cannibalization  cost  was  much  smaller  than  the  other 
costs.  The  reason  for  this  is  understood  by  examining  the  values 
of  the  cannibalization  probability  at  the  end  of  each  time  horizon 
considered.  Table  5.4  contains  the  listings  of  the  cannibalization 

4 

Appendix  IV  contains  the  analytic  form  of  the  multiple 
convolutions  required  for  the  repair/replacement,  and 
cannibalization  cost  calculations. 

^Appendix  V outlines  the  method  of  Gaussian  Quadrature  which 
was  used  to  numerically  integrate  the  repair/replacement  cost, 
cannibalization  cost,  and  the  revenue  function. 


TABLE  5.4  Cannibalization  Probability  as  a Function  of  Time  and 


Spare  Configuration 

I 

1 


T = 
Ov 

5 Years 
m2,  m3) 

D-j  (T) 

d2(t) 

MT> 

( 4 

2 

1) 

2 ?5  X 10‘3 

2.06  X 10"3 

2.93  X 10"6 

( 4 

2 

2) 

9.54  X 10"5 

8.70  X 10‘5 

1.43  X 10"9 

( 4 

3 

0) 

3.19  X 10"2 

9.42  X 10“3 

3.30  X 10"3 

( 5 

2 

0) 

1.58  X 10"2 

5.51  X 10"2 

3.47  X 10"3 

T = 

(mp 

10  Years 
m2,  n^) 

D-|  (T) 

d2(t) 

d3(t) 

( 4 

2 

1) 

1.23  X 10"2 

1.07  X 10"2 

4.33  X 10"6 

( 4 

2 

2) 

9.99  X 10"4 

8.77  X 10"4 

2.91  X 10~9 

( 4 

3 

0) 

1.17  X 10"1 

3.62  X 10"2 

4.86  X 10'3 

( 5 

2 

0) 

6.03  X 10'2 

1.65  X 10"1 

4.39  X 10~3 

T = 
(nip 

15  Years 
m2>  m3) 

D-,(T) 

d2(t) 

°3(t) 

( 4 

2 

1) 

2.72  X 10“2 

2.43  X 10~2 

4.40  X 10"6 

( 4 

2 

2) 

3.23  X 10~3 

2.91  X 10“3 

2.96  X 10“9 

( 4 

3 

0) 

1.83  X 10"1 

6.16  X 10“ 2 

X 

1 

O 

1 

CJ 

( 5 

2 

0) 

9.76  X 10~2 

2.65  X 10_1 

4.43  X 10“3 

r 
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probability,  D.(T),  for  each  spare  configuration  over  the  three 
J 

time  horizons  considered  . 

In  order  to  obtain  the  present  worth  of  the  return  or  revenue 
function,  K^,  one  must  analyze  all  possible  state  transitions 
between  the  "up"  and  "down"  states  for  each  of  the  three  component 
types.  The  analysis  must  include  all  possible  transitions  for  each 
spare  allocation.  For  example,  component  type  1 must  be  analyzed 
assuming  4 or  5 spares,  component  type  2 can  have  2 or  3 spares, 
and  component  type  3 can  have  0,  1,  or  2 spares. 

Because  of  the  assumption  of  statistical  independence  of 
the  respective  failure  and  replacement  times,  stated  in  Chapter  I, 
one  can  obtain  the  unit  availability  function,  P^(t),  by  calculating 
the  respective  component  availability  functions,  i.e.  the  probabil- 
ity of  0,  1 , or  2 components  of  each  type  being  up  at  any  time,  t. 
Define  q..(t)  as  the  probability  of  having  i components 

* J 

of  the  jth  type  being  up  at  any  time  t,  for  i=0,  1,  or  2,  and 
0=1,  2,  and  3.  P^(t)  was  defined  earlier  as  the  time  dependent 

probability  of  having  exactly  i units  up  at  time  t.  Then, 

(5.6)  P^ ( t ) = Pr  (exactly  2 units  are  up  at  time  t.} 

{2  components  of  type  1 are  up,  2 components'] 
of  type  2 are  up,  and  2 of  type  3 are  l 
up  at  time  t.  | 

= qai^)  q22(t)  q23(t). 

6For  each  time  horizon,  T,  ten  values  of  the  cannibalization 
probability  were  calculated;  the  last  being  T x R,n,  where 
R-|q(=0.9869)  was  the  10th  root  used  for  the  Gaussian  Quadrature 
integration. 


S' 
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7 

Similarly  , 
(5.7)  PQ(t) 


Since  there  must  be  0,  1,  or  2 units  up  at  any  time  t,  then, 

(5.8)  P1(t)  = 1 - P0(t)  - P2(t). 

In  order  to  calculate  each  of  the  q..  terms,  a state  transition 

10 

analysis  is  performed  by  modelling  the  transitions  of  spare  con- 
sumptions as  a Markov  Process,  with  an  absorbing  state;  the  ab- 
sorbing state  is  the  state  achieved  after  the  (N+m.-l)th  replace- 

th  ^ 

ment  of  the  j component  type,  i.e.  after  (N-l)  components  have 

been  cannibalized. 


= Pr  {exactly  0 units  are  lip  at  time  t.) 


0 components  of  type  1 are  up,  or 
= Pr  <0  components  of  type  2 are  up,  or_ 

0 components  of  type  3 are  up  at  time  t. 


qqq  +qqq  +qqq  +qqq 

01  02  03  01  02  13  01  12  03  01  12  13 

+ QQQ  +qqq  + q q q +qqq 

Hlr02H03  h11m02h13  m11h12h03  m01H02h23 

+ qqQ  +qqq  + q q q +qqq 

01  22^03  01  22  23  21  01  03  21  02  23 

+ q2iq22q03  + q0iqi2q23  + q0iq2iqi3  + qiiq22q03 

+ qqq  +qqq  +qqq  . 

21  12  03  11  02  23  21  02  13 


^For  notational  convenience,  the  time  dependence  of  q.jj(t) 
will  not  be  shown  in  this,  and  any  subsequent  calculations. 


State  Transition  Analysis 


For  any  given  component  type,  define  S.  . (t)  as  the  time  de- 

1 >K 

pendent  probability  of  being  in  the  state  (i,k);  i.e.  having 

I 

i components  up  with  k spares  remaining.  The  state  transition 
diagrams  for  each  possible  spare  configuration  are  contained  in 
Figures  5.1  through  5.5.'  The  possibility  of  using  a cannibalized 
component  to  replace  a failed  component  is  illustrated  by  the 
final  transition  to  state  (1,-1).  The  negative  index  for  the 
number  of  spares  remaining  indicates  that  the  replacement  came  from 
a cannibalized  component. 

The  last  possible  transition,  to  state  (1,-1),  is  the 
transition  to  a cannibalized  state.  After  cannibalization,  one 
can  not  achieve  full  system  operational  capability  because  the 
very  process  of  cannibalization  eliminates  one  unit  from  active 
service;  consequently,  the  best  capability  one  achieves  after 
cannibalization  is  a working  system  of  1 operational  unit. 

Figure  5.1  contains  the  state  transition  diagram  for  com- 
ponent type  1,  having  5 spares  initially.  There  are  19  possible 
states;  the  state  transition  equations  will  be  transformed  into 
a system  of  19  coupled  first  order  differential  equations. 

Referring  to  Figure  5.1,  the  following  analysis  results: 

I 

During  an  arbitrary  interval  of  duration  dt,  component  type  1 
will  go  from  a state  in  which  2 units  are  up  with  5 spares 
remaining,  state  (2,5),  to  a state  with  1 unit  up  and  5 spares 


Figure  5.1  State  Transition  Diagram  for  a Configuration  with  5 Spares  Initially  Provisioned. 


Configuration  with  4 Spares  Initially  Provisioned. 


Figure  5.3  State  Transition  Diagram  for  a Configuration  with  3 Spares  Initially  Provis 


Figure  5.4  State  Transition  Diagram  for  a Configuration  with  2 Spares  Initially  Provisioned. 


Figure  5.5  State  Transition  Diagram  for  Configurations  having  1 Spare,  and  0 Spares 
Initially  Provisioned. 
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remaining,  state  (1,5),  with  probability  2X  dt.  Thus  the  probabil- 

1 I 

ity  that  it  remains  in  state  (2,5)  is  (1  - 2X,dt).  Once  it  is  in 

1 I ' 

state  (1,5),  it  can  remain  in  that  state  with  probability 
(1  - y-jdt  - Xjdt),  or  it  can  undergo  a replacement  and  go  to 
state  (2,4)  with  probability  p^dt,  or  it  can  fail  with  probability 
Ajdt,  and  go  to  state  (0,5).  The  19  state  transition  equations 
for  component  type  1 follow. 


(5'9)  S2,5(t+dt)  = (1"2Xidt)  S2,5(t) 


S2,5(0)=1-° 


Si>5(t+dt)  = ( 1 -p i d t-X t d t ) Slj5(t)  + 2X-j d t S2>5(t) 
s0,5(t+dt)  » (1-^dt)  SQ>5(t)  + ^dt  S^5(t) 

S2f4(t+dt)  = (l-^dt)  SM(t)  + V,dt  S1>5(.t) 

S]  4(t+dt)  = (1  X dt-p  dt)  S1>4(t)  + y^t  SQ>5(t) 
+ 2X1dt  S2  4(t) 

SQj4(t+dt)  = (1-y^t)  S0>4(t)  + X^t  S1>4(t) 

S2,3(t+dt)  = (1-2X1dt)  S2,3(t)  * V1  S1.4(t) 

S1  3(t+dt)  = (1-X-jdt-y-jdt)  S 1 3(t)  + P^t  S^4(t) 
+ 2X^dt  S2>3(t) 


sQ>3(t+dt)  = (l-yt)  s0>3(t)  + x^t  sli3(t) 


(5.9)  ! 

(cont'd) 


;2  2(t+dt)  = (i-2X1dt)  s2  2(t)  nyt  slf3(t) 

> 2(t+dt)  = (1-X-jdt-y^dt)  S1  >2(t)  + P-j dt  s0#3(t) 

+ 2Xjdt  s2  2(t) 

S0,2(t+dt)  = (1_|ildt)  S0,2(t)  + Xldt  S1 ,2(t) 

S2  1 (t+dt)  = (l-2X1dt)  S2J(t)  +;  y^t  S12(t) 

51  -|(t+dt)  = (l-y^dt-X^dt)  S1  ^t)  + y^t  Sq  ^(t) 

+ 2X^dt  S2>1(t) 

so,i(t+dt)  = (1"yidt)  so,i(t)  +xidt  si»i(t) 

52  Q(t+dt)  = (l-2X1dt)  S2  Q(t)  + y-jdt  S]  #1  (t) 

S1>Q(t+dt)  = (1-^dt)  Slj0(t)  + y-|dt  S0fl(t)  + 2X-j d t S2>0(t) 

so,o(t+dt)  = (1'yidt)  so,o(t)  + Xidt  Sl,0(t) 


S,  . (t+dt)  = y.dt  At) 


’ 


11 


Using  the  definition  of  the  derivative, 

Si,k<t+dt> 


(5.10)  S.  (t)  = lim 
1,k  dt  0 


dt 


the  following  system  of  first  order  differential  equations  may  be 
constructed  for  the  first  component  type  when  it  has  5 spares. 


(5.11)  Sw(t>  = -2  X,S2>5(t) 

s,t5(0  - -U,  + »*i)  s,  5(t)  ♦ 2 X,s2>5(t) 

5O.B(t>  = '“I  S0.5(t)  + XlS1.5(t) 

*2.4(t)  * -2  V2,4(t)  * 


S,  10)  = 1.0 

c9b 


Sl,4(t)  ■ 
*0,4^  = 
^3(t)  = 
^l,3(t)= 
*0.3(t)  = 


-(xn  + Pi)  S,  -(t)  + y,  ft)  + 2 A.S  (t) 


2,4 


1 “V  n,4vw  H1  o,5' 

“yl  S0,4^  + X1  Sl,4^^ 

-2V2,3(t)+Vl,4(t) 


“(X1  +W1)  Sl,3(t)  +yi  S0,4(t)+2X1  $2,3(t) 
^1  S0,3(t)+Xl  $l,3(t> 


$2,2(t)  = 


2A  S (t)  + y.  S (t) 

1 2,2  1 


S,  = “U,  + vO  S (t)  + y S ft)  + 2A  S (t) 


(cont'd) 


S0.2 

(t)  = 

- V 

1S0. 

2(t)  + 

xlsl. 

2(t) 

$2.1 

(t)  = 

-2 ; 

SS2 

.,<t)  + 

*lsl, 

2(t) 

*.,1 

(t)  = 

-(X 

1 + 

Sl.l 

(t)  + 

11  Is 

CO  . 

o 

«• 

(t)  = 

-pi 

so. 

,(t)  t 

Msi. 

l(t) 

S2.0 

(t)  = 

-2  ; 

\..s 
1 2 

ft)  + 

.0 

uisi 

,l(t) 

*1.0 

.(t)  = 

-xi 

si. 

0(t) + 

"1S0, 

,(t) 

*o,c 

,(t)  = 

"yi 

so, 

,o(t)  + 

*1*1. 

,o(t) 

• 

s 

1,- 

,(t)  ■ 

: yi 

so, 

,o‘*> 

From  the  solution  of  these  differential  equations  , it  is 
possible  to  obtain  the  time  dependent  probability  of  0,  1,  or 
2 components  of  the  first  type  being  "up",  or  available,  at 
any  given  time,  q0l , q^,  or  q2i  • 

(5.12)  q0I  = S05(t)  + S0>4(t)  + S0>3(t)  + S02(t)  + S0,(t) 

+ s0>0(t). 


Appendix  VI  contains  the  differential  equations  for  the  other 
configurations,  i.e.  a 2 unit  system  with  4,  3,  2,  1,  or  0 spares 
for  any  component  type. 
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(5.13)  qn  = S1>5(t)  + S1>4(t)  + S1>3(t)  + S1  (t)  + Slfl(t) 

+ SliQ(t)  + S3  (t) 

(5.14)  q21  = S2f5(t)  + S2>4(t)  + S2t3(t)  + Sz2(t)  + S2J(t) 

+ S2.0<‘> 

Once  the  time  dependent  unit  availability  function,  P^(t), 
is  determined,  the  present  worth  of  the  return  or  revenue  function 
is  calculated  using  Gaussian  Quadrature  to  perform  the  following 
integration: 

T 2 

(5.15)  Kd  = / E V,  P- (t)  exp(-rt)  dt, 

* o 1=1  1 1 

for  r = 10%,  15%,  and  20%;  T = 5,  10,  and  15  years;  V-j  = $8000.00 
per  year,  ~ $16000.00  per  year. 

Tables  5.5,  5.6,  and  5.7  contain  the  final  solutions  for  the 
example  problem  for  each  of  the  three  time  horizons,  and  with 
continuous  interest  rates  of  10%,  15%,  and  20%  per  year,  respective- 
ly. In  each  case,  the  allocation  of  4 Line  Printers,  2 Drum 
Plotters,  and  1 Card  Reader  resulted  in  the  minimum  net  present 
worth  of  the  total  costs,  K^.  This  is  primarily  due  to  the  over- 
whelming effect  of  the  respective  procurement  costs,  , in  com- 
parison with  the  other  costs.  In  each  case,  allocation  (4,2,1) 
had  the  minimum  procurement  cost  and  replacement  cost,  while 
allocation  (4,2,2)  had  the  smallest  cannibalization  cost. 

Allocation  (4,3,0)  yielded  the  maximum  revenue  but  it  was  not  large 


TABLE  5.5  Present  Worth  of  the  Total  Costs  for  each  Spare 
Configuration,  with  a 10%  Annual  Interest  Rate. 


. 


K 


0 


= Net  Present  Worth  of 


Total  Operational  Costs  ($). 


T = 

5 Years 
^2  ’^3) 

II 

K1  + 

k2 

+ 

k3  - 

K4 

( 4 

2 

1) 

80843.72 

124550.00 

2451.68 

6.24 

46164.20 

( 4 

3 

0) 

81716.23 

1281 50.00 

2826.15- 

68.63 

49328.55 

( 4 

2 

2) 

86275.25 

130000.00 

2451.79 

0.23 

46176.77 

( 5 

2 

0) 

94466.79 

131100.00 

2979.33 

137.71 

39750.25 

ii  * — • 

E 

1— 

10  Years  KQ 
,m2  ,rn3 ) 

K1  + 

k2 

+ 

k3  - 

K4 

( 4 

2 

1) 

61873.78 

124550.00 

2619.05 

67.55 

65362.82 

( 4 

3 

0) 

62965.42 

128150.00 

3121.49 

563.24 

68869.31 

( 4 

2 

2) 

67246.08 

130000.00 

2620.25 

4.28 

65378.45 

( 5 

2 

0) 

76501.27 

131100.00 

3171.47 

945.56 

58715.76 

T = 
(mr 

15  Years  KQ 
.n^  >1^3) 

K1  + 

k2 

+ 

K3 

K4 

( 4 

2 

1) 

50505.70 

124550.00 

2698.58 

175.62 

76918.50 

( 4 

3 

0) 

52107.65 

128150.00 

3184.96 

1154.69 

80382.00 

( 4 

2 

2) 

55784.68 

130000.00 

2702.56 

15.37 

76933.25 

( 5 

2 

0) 

65924.56 

131100.00 

3236.74 

1894.70 

70306.88 

131100.00  3236.74  1894.70  70306.88 
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TABLE  5.6  Present  Worth  of  the  Total  Costs  for  each  Spare 
Configuration,  with  a 15%  Annual  Interest  Rate. 


Kq  = Net  Present  Worth  of  Total  Operational  Costs  ($). 


T = 5 Years 
(m] ,m2,m3) 

K0 

K1  + 

k2 

+ 

K3  ' 

k4 

( 4 

2 

D 

84883.72 

124550.00 

2197.64 

5.06 

41868.98 

( 4 

3 

0) 

86104.14 

128150.00 

2514.47 

• 

55.71 

44616.04 

( 4 

2 

2) 

90318.05 

130000.00 

2197.73 

0.18 

41879.86 

( 5 

2 

0) 

97772.14 

131100.00 

2649.47 

112.16 

36089.49 

T = 
(mr 

10  Years 

'^2  ,m3 ) 

7^ 

O 

1) 

Ki  + 

K2 

+ 

K3  - 

K4 

( 4 

2 

1) 

71671.83 

124550.00 

2317.32 

46.41 

55241.91 

( 4 

3 

0) 

73032.29 

128150.00 

2728.46 

392.37 

58238.54 

( 4 

2 

2) 

77065.88 

130000.00 

2318.14 

2.88 

55255.14 

( 5 

2 

o) 

85271.24 

131100.00 

2788.60 

663.02 

49280.38 

T = 
(m1 . 

15 
,m2 : 

Years  Kn  = 

,m3) 

Ki  + 

k2 

+ 

*3  - 

K4 

( 4 

2 

1) 

65513.94 

124550.00 

2359.16 

104.06 

61499.28 

( 4 

3 

0) 

67158.43 

128150.00 

2761.93 

710.16 

64463.66 

( 4 

2 

2) 

70858.67 

130000.00 

2361.46 

8.75 

61511.54 

( 5 

2 

0) 

79560.11 

131100.00 

2825.07 

1172.81 

55537.77 

t 
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TABLE  5.7  Present  Worth  of  the  Total  Costs  for  each  Spare 

I 

Configuration,  with  a 20%  Annual  Interest  Rate. 


Kg  = Net  Present  Worth  of  Total  Operational  Costs  ($). 


T = 
(lH], 

5 Years 
m2»m3) 

ii 

o 

K1  + 

K2 

+ 

K3  " 

K4 

( 4 

2 

1) 

88377.99 

124550.00 

1975.37 

4.11 

38151.49 

( 4 

3 

0) 

89897.16 

1281 50.00 

2244.02 

45.27 

40542.13 

( 4 

2 

2) 

93814.66 

130000.00 

1975.44 

0.15 

38160.93 

( 5 

2 

0)  1 00630. A6 

131100.00 

2362.97 

91.46 

32923.97 

T = 10  Years 
(m1 ,m2,m3) 

ii 

o 

K1  + 

k2 

+ 

K3  - 

k4 

( 4 

2 

1) 

79130.12 

124550.00 

2061.31 

32.14 

47513.33 

( 4 

3 

0) 

80739.18 

128150.00 

2399.60 

275.54 

50085.96 

( 4 

2 

2) 

84539.28 

130000.00 

2061 .87 

1.96 

47524.55 

( 5 

2 

0) 

91890.06 

131100.00 

2464.05 

468.94 

/ . 

42142.93 

T = 
(r^ 

15  Years 
1^2  >^2  ^ 

II 

o 

K1  + 

K2 

+ 

K3 

f 

h 

( 4 

2 

1) 

75780.88 

124550.00 

2082.81 

63.04 

50914.97 

( 4 

3 

0) 

77580.22 

128150.00 

2416.77 

477.18 

53463.73 

( 4 

2 

2) 

81163.89 

130000.00 

I 

2084.16 

5.08 

50925.35 

( 5 

2 

0) 

88804.84 

1 31100.00 

2484.90 

744.15 

45524.21 
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enough,  in  comparison  with  the  other  costs,  to  alter  the  optimal 
spare  configuration. 

To  see  the  effect  of  the  stochastic  nature  of  the  revenue 
or  return  function,  one  needs  only  to  calculate  upper  and  lower 
bounds  on  the  expected  revenue  by  assuming  P^(t)  = 1,  for  all 
time.  Then,  a lower  bound  on  the  revenue  function  can  be  cal- 
culated by  assuming  1 unit  was  operational  throughout  the  entire 
time  horizon  yielding  an  annual  revenue  of  $8000.00  per  year; 
likewise,  an  upper  bound  can  be  calculated  by  assuming  both 
units  were  operational  throughout  the  entire  time  duration.  With 
the  relatively  small  MTTR  values,  in  comparison  with  the  respective 
MTBF  values,  one  would  expect  the  present  worth  of  the  revenue 
function  to  lie  between  the  bounds  K^,  and  K^U , w here  the 
respective  bounds  are  calculated  as  follows: 


(5.16)  K4L  = 


(5.17)  K4U  = 


/V  exp  (-rt)  dt 
o 1 

V-j  (1  - exp(-rT))/r 
$8000  (1  - exp(-rT) )/r. 

exp  (-rt)  dt 

o 2 

V2  (1  - exp(-rT))/r 
$16000  (1  - exp(-rT) )/r. 


Table  5.8  contains  values  of  the  upper  and  lower  bounds  on 
the  present  worth  of  the  expected  revenue  function  for  interest 
rates  of  10%,  15%,  and  20%  per  year,  over  time  horizons  of  5,  10, 


r 
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quirements.  It  is  a present  worth  analysis,  so  the  derived  cost 
■ comparisons  are  only  valid  if  compared  over  equal  time  horizons. 

An  example  problem  was  presented  in  which  there  were  two 
active  constraints:  a performance  constraint,  which  specified  the 
minimum  required  probability  of  spare  adequacy,  and  a resource  con- 
straint, which  imposed  a limitation  on  the  total  investment  in 
spares.  The  complexity  of  the  cost  calculations  was  formidable 
for  all  but  the  smallest  system  configurations;  this  is  mainly  due 
to  the  complexity  of  the  expression  for  the  cannibalization  prob- 
ability, products  of  (P-1)  double  integrals  for  each  of  the  P com- 

y j 

ponent  types.  The  problem  is  further  complicated  by  the  necessity 

l 

of  calculating  the  multiple  convolutions  for  the  marginal  density 


[ 


functions  of  the  terminal  failure  times. 

Gaussian  Quadrature  proved  to  be  an  efficient  numerical 

I 

technique  for  performing^ both  the  multiple  integrations  for  the 


• I 


time  dependent  cannibalization  probability,  and  the  single  integra- 


i I 

>.  i 


tions  for  each  of  the  respective  present  worth  calculations. 


CHAPTER  VI 


CONCLUSIONS  AND  SUGGESTIONS  FOR  FURTHER  RESEARCH 


The  primary  objective  of  this  research  was  to  develop  an 
operational  cost  model  which  allowed  one  to  analyze  the  effects  of 
cannibalization  on  spare  inventory  allocations.  The  model  was  not 
restricted  to  areas  in  which  instantaneous  replenishment  was  assumed; 
it  was  developed  under  the  assumption  of  exponential  failures  and 
a general  replacement  density  function. 

A mathematical  formulation  of  the  model  is  restated  below: 
Minimize  K0  = ^1  + K2  + K3  " where  1s  the  Present  wortfl  of 
the  net  operational  costs,  K-j  is  the  present  worth  of  the  manufac- 
turing or  procurement  costs,  l<£  is  the  present  worth  of  the  repair/ 
replacement  costs,  K-j  is  the  present  worth  of  the  cannibalization 
costs,  and  is  the  present  worth  of  the  return  or  revenue  function. 
The  equations  for  each  of  the  major  terms  are  defined  as  follows: 

P P 

(6.1)  K.  = N £ a-  + £ a-m.. 

1 j=l  J j=l  J J , 

T P N+m  -1  /in  m \ 

(6.2)  K = / £ b . (t)  £J  i(V,;*g.(t))  * R.U'(t)  exp(-rt)dt. 

2 o j=l  J i=l  J J J 

T P 

K = / (N-l)  £ c.D,(t)  exp(-rt)dt,  where 

3 o j=l  J J 


(6.3) 
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D.(t)  is  the  time  dependent  cannibalization  probability,  defined  as 

J 

follows: 


(6.4)  D.(t)  = n Si  1}(x)  h ^+nlk  ^(y)  dx  dy,  where 

J k=l  o o J K 

m 

. (N-Hnj  1S  the  marginal  density  function  of  the  time  to  the 

J 

xi  XL 

(N+m.-l)tn  failure  of  the  j component  type;  the  analogous  defini- 
J ( 1 ) 

tion  holds  for  h^  k '(t).  The  revenue  function  is  expressed  as 
follows: 


T N 

(6.5)  K . = / Z V.  P.(t)  exp(-rt)dt. 
q o i=o  1 ‘ 


The  general  model  is  also  applicable  when  there  are  one  or  more 
constraints,  which  usually  provide  bounds  on  the  permissible  spare 
allocations.  A constraint  on  the  minimum  required  spare  adequacy 
results  in  a lower  bound  on  the  number  of  spares,  while  an  invest- 
ment constraint  usually  provides  an  upper  bound.  Other  applicable 
constraints  include  a constraint  on  the  total  weight  or  volume  of 
a specific  spare  allocation,  a requirement  on  the  system's  interval 
availability,  and  a maintainability  constraint  on  system  downtime. 

A specific  example  problem  was  solved  which  illustrated  the 
applicable  cost  calculations.  Gaussian  Quadrature  was  used  to 
perform  the  required  integrations. 

Because  of  the  complexity  of  the  marginal  density  functions  of 
the  terminal  failure  times  for  large  system  configurations,  a general 
allocation  model  was  not  computationally  feasible;  however,  one 


should  be  able  to  model  his  allocation  problem  using  the  methodology 
developed,  and  proceed  to  a complete  solution  using  the  same  steps 
presented  in  the  example  problem. 

The  methodology  developed  in  this  research  can  be  put  to  use 
in  either  of  two  contexts.  First,  it  provides  a tool,  based  on 
equipment  reliability  and  maintainability  data,  for  an  accurate 
prediction  of  net  operational  costs  on  a present  worth  basis. 

Second,  the  methodolcgy  provides  a mechanism  for  showing  the  re- 
lationship between  a spare  stockage  level  and  the  corresponding 
risk  of  a critical  spare  shortage.  Consequently,  this  analysis 
makes  it  possible  to  study  the  interaction  between  net  costs  and 
the  possible  shortage  risks  over  the  system's  lifetime. 

The  application  of  this  cost  model  can  benefit  both  industrial 
and  military  procurement  strategies.  The  analysis  is  especially 
useful  when  advancing  technology  hastens  the  process  of  obsolescence 
of  current  equipment.  Management  can  utilize  this  model  to 
determine  the  best  of  several  alternative  spare  allocations, 
based  upon  the  present  worth  of  the  net  operational  costs,  assuming 
the  interest  rate  and  time  horizon  are  known.  Consequently,  short 
run  decisions  can  be  made  regarding  the  effectiveness  of  purchasing 
a new  system,  versus  keeping  the  old.  The  present  worth  analysis 

allows  management  to  make  such  a decision  if  the  cost  comparison 

I 

is  made  over  the  same  time  ^interval. 

The  model  is  also  useful  in  situations  having  stringent 
requirements  on  spare  adequacy,  i.e.  aQ  very  close  to  1.0.  Spare 
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provisioning  for  military  weapon  systems  falls  into  such  a category. 

As  in  most  research  endeavors,  there  is  always  room  for 
additional  work;  this  research  is  no  exception.  Some  further 
extensions  of  this  model  are  outlined  below. 

There  are  many  situations  in  which  the  continuous  annual 
interest  rate  is  not  deterministic,  i.e.  it  may  be  a random  variable 
over  a specific  time  horizon.  This  is  especially  true  during  a 
volatile  money  market,  when  the  Federal  Reserve  Board  may,  on  a 
relatively  short  notice,  increase  the  going  interest  rate  charged 
its  member  banks.  Such  an  action  has  a "domino"  effect  in  the 
financial  world,  with  banks  subsequently  raising  the  prime  rate 
charged  their  most  credit-worthy  customers.  Consequently,  a 
stochastic  analysis  of  the  present  worth  of  the  net  operational 
costs  under  cannibalization,  which  considers  the  random  nature  of 
the  interest  rate,  would  be  a valuable  contribution  to  the  state 
of  the  art.  One  may  be  able  to  determine  a probability  density 
function  of  the  present  worth  of  the  net  operational  costs,  and 
develop  risk  profiles  which  illustrate  the  probability  of  the 
present  worth  exceeding  a specified  level.  This  analysis  would 
be  extremely  beneficial  if  applied  prior  to  the  procurement  of  a 
complete  system  or  spare  configuration. 

Another  extension  of  this  model  would  be  to  consider  the  effects 
of  taxes  and  inflation  on  the  optimal  spare  allocation,  under  the 
same  assumptions  used  herein.  With  this  analysis,  the  model  could 
be  tailored  to  the  specific  industries  which  are  subject  to  known 
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tax  rates.  Inflation  is  another  aspect.  Without  considering 
inflation,  a company  may  experience  a cash  shortfall  equal  to  the 
amount  by  which  an  equipment's  accumulated  depreciation  fails  to 
cover  its  replacement  cost.  This  possibility  should  be  considered 
when  procuring  a new  system  or  its  spare  configuration;  a model 
which  includes  the  inflationary  effects  on  procurement  and  replace- 
ment costs  would  be  valuable. 

A further  extension  of  this  basic  model  is  a consideration 
of  catastrophic  failures,  i.e.  where  the  failure  of  one  or  more 
components  causes  the  failure  of  the  entire  unit.  This  was  not 
considered  in  the  analysis  contained  herein,  because  of  the 
assumption  of  statistically  independent  failure  times.  A cata- 
strophic failure  analysis  is  applicable  in  hostile  military 
environments,  where  enemy  actions  directed  against  a single  com- 
ponent or  subsystem  destroy  the  complete  unit.  Consequently, 
cannibalization  of  that  unit  is  prevented. 

Continuing  in  the  same  area,  one  may  be  able  to  develop  a 
cannibalization  cost  model  when  all  the  components  are  subject 
to  multiple  types  of  failure  or  when  a unit  is  permitted  to  operate 
under  a degraded  level  of  performance.  Such  an  analysis  necessitates 
a modified  failure  density  function  for  each  component  type. 

The  present  worth  analysis  of  cannibalization  in  determining 
spare  provisioning  requirements  has  been  an  enlightening  area  of 
research.  The  author  hopes  that  this  work  will  help  motivate 
others  to  continue  in  the  same  area. 
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APPENDIX  I 


CALCULATION  OF  THE  DENSITY  FUNCTION  OF  THE  ORDERED  FAILURE  TIMES 


(i).  , th 

The  term  f . (t)  is  the  density  function  of  the  i ordered 

J th 

failure  time  of  the  j component  type;  it  is  obtained  using  the 
methods  of  order  statistics.  In  general,  the  probability  density 
function  of  the  q**1  order  statistic  from  a random  sample  of  size  L, 
from  a continuous  population,  x,  -°°  < x < ® is: 


(Al.l)  p(q)(x)  = 


q-1.  1,  L-q 


(P(x))q_1  p(x)  (l-P(x))L~q 


(q-1)!  1!  (L-q)! 


(P(x))q_1  p(x)  (l-P(x))L"q 


where  p(x)  is  the  continuous  probability  density  function  of  the 
random  variable  X,  and  P(x)  is  the  corresponding  cumulative 
distribution  function: 


(A1.2)  P(x)  = Pr  (X  < x). 


For  the  general  system  configuration  of  Figure  3.1  (p.  21), 

f.^(t)  is  calculated  belov/: 

J 


r 


1 


y 


A 
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(A1.3)  fj(l)(t) 


N 


^i-1 , 1,  N-i  l 3 3 

for  i = 1 , 2,  ...»  N,  and 

.1-1 


(F,(t))1_1  f,(t)  (l-F-(t))N_1 
! J 


(1)  / N4mi  i i-i  i • 

f'  (t)  = J (F.(t))1  f.(t)  (T -F . ( t ) ) J , 

3 \1-1,  1,  N+m.-i/  3 3 ' J 

J 


N+m.i-1 


for1  i = N+l , N+2 N+m.-l. 

J 


The  term  f.(t)  is  the  probability  density  function  of  the  unordered 
J th 

time  to  failure  of  the  j component  type;  F.(t)  is  the  cumulative 

J 


probability  distribution  function  of  the  unordered  time  to  failure 


.th 


of  the  j component  type. 


(1), 


(1) 


Using  Equation  (A1.3),  the  expression  for  f.v  '(t),  as  well  as 

J 


R.  (t),  is  derived  as  follows: 
J 


fll  N o N-l 

f.U,(t)  = F (t)  f (t)  (l-F.(t)) 

3 \ 0 , 1,  N-l  / 3 3 3 


(A1.4) 


= Nf j (t)  (l-Fj(t)) 


N-l 


= NX.  exp(-Xjt)  (exp(-Xjt)) 


N-l 


= NXi  exp(-NX.t) 

J J 


= 3.  exp(-B-t),  where  3.  = NX.. 
J J J J 


(A1.5) 


(1) 


R.'  ' (t)  = / f.(1)(t)  dt 
3 t 3 


= exp(-3jt) 


^The  expression  for  ,f  -^(t),  when  i > N and  f.-(t)  is  not 
exponential  density  function,  is  an  approximation.  It  is  based 


an 

based  on 

the  fact  that  the  maximum  sample  size  is  N+m.,  even  though  no 
more  than  N components  are  operating  at  any  iime.  For  the  expo- 
nential it  is  exact. 


J 


r 


APPENDIX  H 


TABLES  OF  THE  NORMALIZED  PRESENT  WORTH  OF  THE  REPAIR/REPLACEMENT  COSTS 
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(A3. 9)  f ' '(t)  « 2X,(1  - e * ) e 


X3l\  "1X3t 


Substituting  the  correct  forms  of  Equations  (A3.7-A3.9)  into 
Equations  (A3.4-A3.6),  and  setting  the  upper  limit  of  the 
integrations  to  infinity,  yields  the  following: 


(A3. 10)  D («)  = I / 24\  A (1-e  XlX  )2  e 2XlX  (1-e  X^)  e ^ dxdy 

1 o o 1 2 

oo  y -a-iX  2 “2^-jx  "A^y  x3y 

X/  .C  24A  A (1-e  ) e (1-e  ) e dxdy 

0 0 13 

(A3. 11)  D2H  = £ )jy24A1A2(l-e  ^ ) e U2X  (1-e  X]y)  e 2XlY  dxdy 


(A3. 12)  D3(~) 


• y ..  . , -a9x  -ia,x  -A,y.  -i'A,y 
X / / 4Xo*o(l-e  2 ) e 2 (1-e  3 ) e 3 dxdy 


oo  v -A,x  -lAoX  -A.y  ^ ~2A.y 

= / r 24A,A0(l-e  3 ) e 3 (1-e  1 ) e 1 


oo  13 


X I*  / 4A2A3(l-e'X3X  ) e"1X3  (l-e'X2y)  e"U2y  dxdy 
0 0 £ ^ 


LI.  < 
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Performing  the  required  double  integrations,  and  multiplying, 
yields  the  steady  state  cannibalization  probability  as  a function 
of  the  respective  hazard  rates,  X^,  X^,  and  X3» 

(A3. 13)  (°°)  = {6X2(1/2X2  - 6/(2X1+lX2)  + 6/(2X1+2X£)  - 3/(4X1+lX2) 

+ 3/(4X1+2X2)  - 8/(3X1+2X2)  + 8/ ( 3X1 +1 X£) ) > 

X { 2X3 (1/2 X3  - 6/(2X1+lX3)  + 6/(2X1+2X3)  - 3/(4X1+lX3) 

+ 3/(4X1+2X3)  - 8/(3X1+2X3)  + 8/ (3X] +1 X3 ) ) } 

(A3. 14)  D2(«)  = (12X1(1/12X1  - 2/(2X1+X2)  + 1/(2X]+2X2)  - 2/(3X1+2X2) 

+ 2/(3X1+lX2)  - 2/(4X1+lX2)  + 1 / (4X] +2X£) ) } 

X{2X3( I/2X3  ' 2/(lX2+n3)  + 2/(lX2+2X3) 

- 1/(2X2+2X3)  + 1/(2X2+1X3))} 

(A3. 15)  D3(«o)  = { 1 2X-j  (1/12X-J  - 2/(2X1+X3)  + 1/(2X1+2X3)  - 2/(3X-,+2X3) 

+ 2/(3X1+X3)  - 2/(3X1+2X3)  + 1/(4X1+2X3) )) 
X{2X2(1/2X2  - 2/(X2+X3)  + 2/(2X2+X3) 

- 1/(2X2+  2X3)  + 1/(X2+2X3))} 
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APPENDIX  IV 

CALCULATION  OF  MARGINAL  DENSITY  FUNCTIONS  OF  TERMINAL  FAILURE  TIMES 

. 

In  the  example  problem  of  Chapter  V,  the  following  probability 

density  functions  are  implicitly  used  in  the  calculation  of  the 

cannibalization  probability:  h^^(t),  h^6^(t),  h2^(t),  h2^(t), 

h3(,)(t),  h3<2>(t),  and  h3(3)(t). 

The  term  h^^(t)  is  the  probability  density  function  of  the 

time  to  the  fifth  failure  of  the  first  component  type.  Because  the 

replacement  time  is  non-zero,  the  required  density  function  is  the 

density  of  the  sum  of  the  sum  of  5 ordered  failure  times,  and  4 

replacement  times. 

tc)  fll  *5  *4 

(A4.1)  h/&V)  = fl  ;(t)  *g,(t)  , 

ii  ' i 

(1)  *5 

where  f (t)  is  the  fifth  multiple  convolution  of  the  probability 

density  function  of  the  first  ordered  failure  time  of  the  first 

*4 

component  type;  g-j(t)  is  the  fourth  multiple  convolution  of  the 
first  component's  replacement  density  function. 

The  term  hj^(t)  is  defined  as  the  probability  density  function 
of  the  time  to  the  sixth  failure  of  the  first  component  type.  It  is 

I (5) 

calculated  in  a way  similar  to  h^  v /(t). 

(6)  fll  *6  *5 

(A4.2)  h]  b (t)  = fj  (t)  *y](t)  . 


<*-**:■ 


nuili.iiyii.lJ  , i ip  . iiiipiipwww^ippipwiwp^ipi  ■•  ■-  i— ' — v ■«.  wfii  > n 
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(3) 

The  term  (t)  is  the  probability  density  function  of  the 
time  to  the  third  failure  of  the  second  component  type. 

(3)  (1)  *3  *2 

(A4.3)  h^  V)  = V (t)  * 92(t)  ’ 

(4) 

The  term  h^  ( t)  is  the  probability  density  function  of  the 

time  to  the  fourth  failure  of  the  second  component  type. 

(4)  (1)  *4  *3 

(A4.4)  h2  (t)  = f2  (t)  * g2(t)  . 

The  term  h^^(t)  is  the  probability  density  function  of  the 
time  to  the  first  failure  of  the  third  component  type.  It  is  cal- 
culated, using  the  method  of  order  statistics,  as  follows: 

(A4.5)  >3(1,(t)  = rj  F3(t)°  f3(t)  (1-F3(t)) 

- 2 f3(t)  ()-F3(t)) 

= 2 X3  exp(-X3t)  (l-(l-exp(-X3t))) 
s 2X3  exp(-2X3t), 

(2)  (3) 

The  terms  h,  (t),  and  h , are  the  density  functions  of 

O O 

the  second  and  third  failures,  respectively,  of  the  third  component 
type. 

(A4.6)  h3(2)(t)  = f3(1)(t)  2 * g3(t). 

Ml  (1)  *3  *2 

(A4.7)  h3'  V)  = Y (t)  * 93(t)  • 


I 


! 
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Under  exponential  failure  and  replacement  density  functions,  the 
respective  marginal  density  functions  of  the  terminal  failure  times 
are  calculated  as  follows;  the  calculations  rre  based  on  the  defini- 
tion of  the  convolution  operation  given  in  Equation  (3.5)  (p.  26). 


(A4.8) 


= y '(t)  *9,(t) 


<'V*5 

^ t .%3 


1 


= B^Bjt)  expj-^t)  * p^pjt)  expt-p-jt) 


3! 


- eA.4/(4!  3!)  < 6t  exp(-B-t) 

l 1 

• 3 3 

- 106t  exp^t)  - 24t  exp(-p  t) 

a 5 a 5 * 

2 1 2 1 
+ 720t  exp(-B  t)  - 360t  expt-p.t) 

~V  °16 

- 2880t  exp(-B  t)  - 2160t  exp(-p.t) 

“l7  V . 

+ 5040  expt-Bjt)  - 5040  expf-p^)  I , 


a.|8 


a^8 


where  = p^  - Bj » and  Bj  = NX^ 


J/ 
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(A4.9)  h/^U)  - f/^U)  *9,(0 

= P^^t)5  expjJS.il)  * (p1  t)4  expC-v^t) 

= B.V5  /(5!  4!)  | 24  t5  exp(-g  t) 

\ -s»-^ 

- 600  t exp (-8  t)  +120  t exp^t) 

~V  1 °1 

+ 7200  t3  exp(-Bi  t)  + 2830  t3  expj-y^) 

2 ^ 2^ 

- 50400  t ex^(-g-|t)  + 30240  t exgl-j^t) 

• “l  “l 

+ 201600  t exp(-git)  + 161280  t expf-p^) 
°1  °1  ^ 

- 362880  ex^-Bit)  + 362880  ex^-y^t) 


o a, 

1 1 


1 


where  ct^  = y-j  - 3-j » and  = Nxv 


(A4.10)  h?(3\t) 


f20)(t)*3  * 92(t) 


*2 


= B2(B2tr  exp(-y2t) 

= 8 3 y?2/(2!)  It2  ex^B^  " 4t  1*^2^ 

A q22  “2 

-2t  expt-V-t)  + 6 


“23  1 

6 exp(-P?  t ' )y , 

a.4  J 


“24 


where  a2  = p?  - B2,  and  B2  = NX^ 


w - 


i.  • 


(A4.ll) 


h2(4,(t, 


- f <’ V4  * 


2 (t)  * g2(t) 


*3 


&2(B2t)  §X£izl2i)  * P2(p2 ** 

3 i _ 21 


= ^2  v2  ^ ^2  exp(- 


S<=M> 


-18  exp(-32t)  + 6 exp(-y2tl 

°2  V 


+72  t exp^-Bot)  + 48  t exg(-y^t) 


- 120  exp(^g^t)  + 120  exp^-y^t) 


where  a„  = P2  - 82>  an(*  82  = NX^. 


(A4.12) 


h3(2,(t) 


0)  *2 

V (t)  * g3(t) 


= B (03t)  exp(-03t)  * y3  exp(-y3t) 
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&3  P3  Jt  exp(-B3t)  + e_xp(-y3t)  - exp(-B3t)l 


a 


where  a = y3  - B3>  and  33  = NX^. 


J 


APPENDIX  V 


NUMERICAL  INTEGRATION  USING  GAUSSIAN  QUADRATURE 

A detailed  analysis  of  Gaussian  Quadrature  may  be  found  in 
Carnahan,  et  al.  [8]  . However,  an  outline  follows,  in  which  the 
details  are  illustrated  for  the  double  integrations  necessary  to 
calculate  D^(t). 

(A5.1)  D,(t)  = d]2(t)  d?3(t), 

where,  assuming  (m^,m2,m3)  = (5,2,0), 

(A5.2)  d12(t)  = / / h/6)(x)  h2^3)(y)  dx  dy,  and 

t v (6)  (1) 

dl3(t)  = £ r h]  (x)  h3  (y)  dx  dy. 

Essentially,  Gauss-Legendre  Quadrature  involves  converting 

integrals  of  the  form  /V(x)  dx,  into  the  form  k /V(z)  dz,  by 

d 0 

a suitable  change  of  variables;  the  latter  integral  is  ultimately 

M 

transformed  into  the  form  of  a finite  summation,  E U.  f(z;). 

i=l  1 1 

Specifically,  let  z=(x-a)/(b-a),  then  x=(b-a)z  + a,  and  dx=(b-a)dz; 
thus, 

(A5.3)  fb f(x)  dx  = (b-a)  f((b-a)z+a)  dz. 

a o 

Next,  the  integral,  £ 'f ((b-a)z+a)  dz  is  transformed  into  a finite 
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summation, 

1 M 

(A5.4)  r f ( (b-a)z+a)  = E U.  f ( (b-a)R.  + a), 

o i=1  1 1 

where  the  R^  terms  are  the  roots  of  the  order  Legendre  Polynomial 
on  the  interval  [0,1]  , and  U.  is  the  corresponding  weighting 
factor  for  the  i root.  The  roots  and  weights  of  the  m"  order 
Legendre  Polynomial  are  given  below. 


Index 

Root 

Weight 

i 

Pi 

Ui 

1 

0.0130467357 

0.0333356721 

2 

0.0674683166 

0.0747256745 

3 

0.1602952158 

0.1095431812 

4 

0.2833023029 

0.1346333596 

5 

0.4255628305 

0.1477621132 

6 

0. 5744371694 

0.1477621132 

7 

0.7166976970 

0.1346333596 

8 

0.8397047841 

0.1095431812 

9 

0.9325316833 

0.0747256745 

10 

0.9869532642 

0.0333356721 

To  perform  the  double  integration  for  d^2(t),  using  Quadrature, 
one  makes  the  following  transformation:  (M=10,  for  a 10^  order  fit) 


M 


(3)  M ffil 

(A5.5)  d (t)  = t E U;(R.t)-h  v '(R.t)  E U .- h/  7 (R- (R.t)  ). 
U i=l  1 1 d 1 j=l  J 1 J 1 


APPENDIX  VI 


STATE  TRANSITION  ANALYSIS  FOR  CONFIGURATIONS  WITH 
4,  3,  2,  1,  OR  0 SPARES 


From  Figure  5.2  (p.  104),  the  following  differential  equations 
may  be  constructed  for  a configuration  with  4 spares. 


(A6.1 ) 


Sl,4^ 


SQ.A™ 


1,2' 


S (t) 

J.l 

so,i 


2,0 


M.o 


0,0 


-2XS2.4(t) 

■(x+p)s1>4< 

— , 4 ( t ) + X$1  4(t) 


s2>4(0)  = 1.0 


-(X+p)SM(t)  + 2AS2>4(t) 


S2.3(t) 

V3(t) 


-2XS2.3(t)+US1.4(t) 


(A+y)S  (t)  + yS  (t)  + 2AS  (t) 
1,3  0,4  2,3 


So.s^) 

>> 

'(t) 


-yS0,3(t)  + XS1>3(t) 


S2.2(t) 


= -2XS  „(t)  + uS,  ,(t) 


2,2 


1.3 


-(X-m)S  (t)  * IIS  (t)  + 2XS  (t) 

1,2  0,3  2,2 


S0,2(t) 


-*jS0.2tt)  + XS1.2(t) 


S2,l(t) 


-(A+y)S  (t)  + yS  (t)  + 2AS  (t) 

1,1  0,2  2,1  ' 


(t)  = 

-ySo, 

,(t)  + 

xsi,i(t) 

(t)  = 

-2\S2 

'.o'1’  ' 

tMSi,i(t) 

(t)  = 

"XS1, 

.0<‘>  + 

ys0>1(t)  + 2xs 

(t)  = 

■ySo, 

0<‘>  + 

xs),0!t> 

(t)  = 

: ^Sn 
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The  time  dependent  probability  of  having  0,  1,  or  2 components  of 
the  first  type  being  up  at  any  time  is  calculated  as  follows: 

(A6-2)  ,0,  = S0>4(t)  + S0  3(t)  ♦ S0  2(t)  ♦ S0>,(t)  + S0  0(t). 

(A6.3)  qjj  = s1>4(t)  + Si  ,3^)  + Sl,2^t^  + Sl,l^  + S1 ,0^ 

+ 5,^(0. 

(A6'4)  q21  = S2,4(t)  + S2,3(t)  + $2,2(t)  + S£,l(t)  + S2,0(t)‘ 

From  Figure  5.3  (p.  105),  the  following  differential 
equations  may  be  constructed  for  a configuration  having  3 spares. 


(A6.5) 


- 

-2XS2>3(t) 

S)3(t)  = 

-(X+y)S!  3(t:)  + 

2AS2.: 

y(t)  - 

-vS0>3(t)  + XSlt 

3<‘) 

S2,2(t)  ■= 

-2XS2>2(t)  + yS1 

.3«> 

V‘>  ‘ 

-(X+y)S1  2(t)  + 

"S0,3 

-ySQ^tt)  + xsi. 

2(t) 

II 

4J 

m 
CM 
• C/> 

-2XS2>1(t)  + yS1 

>2(t) 

Sij(t)  = 

— ( X+y ) S i (t)  + 

vS0.2 

^0,1 = 

-ySo,i(t)  + xsi. 

,(t) 

s2>0(t)  = 

-2XS2>0(t)  ♦ PS, 

,(t> 
* 1 

*l,0(t)  = 

-^l.oW  + ySo, 

,(0 

so,o^  = 

-yS0>0(t)  + ASlf 

0<t> 

Sl.-l(t>  ■ 

: ySo,o(t) 

S2,3<°> 


= 1.0 


>(t) 


r 


i 


If  component  type  2 has  a configuration  with  3 spares,  then  its 
appropriate  state  differential  equations  are  given  in  (A6.5),  and 
the  following  equations  describe  the  time  dependent  probability 
of  0,  1,  or  2 components  of  the  second  type  being  up  at  any  given 
time. 

(A6.6)  qQ2  = Sg^Ct)  + ^0,2^^  * ^0,1^^  + ^0,0^^’ 

(A6.7)  q1?  = S1>3(t)  t Su2(t)  + Slfl(t)  + Q(t)  + S^  ^t). 

(A6.8)  q22  = S2  3(1:)  + S2>2(t)  + S2j(t)  + S2>Q(t). 

Figure  5.4  (p.  106)  contains  the  state  transition  diagram 
for  a configuration  with  2 spares,  which  results  in  the  following 
system  of  differential  equations. 


(A6.9)  S2>2(t)  = "2XS?  2(t) 


S2>2(0)  =1.0 


S1>2(t)  = -(X+u)S1>2(t)  + 2AS2f2(t) 

%2(t)  = "yS0,2(t)  + XSl,2(t) 

S2,i  (t)  = -2AS?  J(t)  + mS]  >2(t) 

Si,i(t)  = -(A+iOS^U)  + yS0>2(t)  + 2AS2J(t) 

‘ ■,jS0.1(t>  * XSl.l(t) 

S2,0^  “ "2XS2,0^  + 

S)>0(t>  * + ^0>)(t)  + 2XS2j0{t) 

S0,0^  ' llSf',0^t^  * “l.O^ 

Vl(t)  * "so,o(t) 


I 
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t 

If  component  type  2 has  a configuration  with  2 spares,  then  the 
following  equations  describe  the  time  dependent  probability  of 
0,  1,  or  2 components  of  the  second  type  being  up  at  any  given 
time. 

(A6.10)  q Q2  = + ^o,i(t)  + S0,0^' 

(A6.ll)  q12  = S1>2(t)  + Sl  ^t)  + S1>0(t)  + S1  ^(t). 

(A6.12)  q^  = S2>2(t)  + S2J(t)  + S2>0(t). 

Similarly,  when  component  type  3 assumes  a configuration 
with  2 spares, 

(A6.13)  q03  = S02(t)  + s0>1(t)  + S00(t). 

(A6.14)  qu  = S]>2(t)  + Slsl(t)  + SKQ(t)  + S^^t) 

(A6.15)  q23  = S2>2(t)  $2  ^2  0^^" 

(One  need  only  use  A3,  and  y3,  in  place  of  A2,  and  pg*  in  t*ie 

solution  of  the  S.  (t)  equations) 

1 , K 

Figure  5.5  (p.  107)  contains  the  state  transition  diagram 
for  configurations  having  1 spare  and  0 spares  initially 
provisioned.  When  component  type  3 has  1 spare,  the  following 
differential  equations  apply. 


3 

t 
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(A6.16) 

s ,(t)  = 
.2,1V 

-2XS2J(t) 

oo 

ro 

«• 

o 

ii 

• 

o 

slj(t)  = 

-(X+y) S1 #1 (t)  + 

2XS2>1(t) 

^0,1^  = 

-^o.i (t)  + xsi. 

,(0 

*2,0^  = 

-2XS2>0(t)  + yS1 

,iW 

^1,0^)  = 

-XS1>0(t)  + ySQj(t)  + 2XS2>Q(t) 

^0,0^  = 

-^S0,0(t)  + XS1, 

0(t) 

^1,-1 (t)  = 

so,o(t) 

(A6.17) 

q03  = S0,l 

(t)  + s0f0(t) 

(A6.18) 

q13  = Sl,l 

(t)  + S1>0(t)  + 

Sl.-ltt) 

(A6.19) 

q23  = S2,l 

(t)  + S2>0(t) 

With 

0 spares. 

the  following  equations  apply. 

(A6.20) 

s2>0(t)  = 

-2*S2j0(t) 

S2,0<°>  ■ 1 

Sl,0(t)  ■ 

-«i,o<t>  + 2XS; 

*0,0^  " 

-^o,o(t)  + xsi 

,0^) 

V,<‘> ! 

! vWl) 

(A6.21 ) 

,03  = S0,( 

3(t) 

(A6.22) 


q13  = Sl,0(t)  + Si,-1 


(t) 


(A6.23) 


q23  S2,0^ 


146 


APPENDIX  VII 

COMPUTER  PROGRAM  FOR  EXAMPLE  PROBLEM  OF  CHAPTER  V 

1.  First,  it  reads  MM,  and  NN,  the  respective  orders  of  the 
integration,  where  MM  is  the  order  of  the  outer  summations, 
and  NN  is  the  order  of  the  inner  summations.  Format(I2,lX,I2) . 

I 

I 

2.  Next,  it  reads  MM  cards  each  with  the  Roots  and  Weights  of 
the  Quadrature  factors  in  the  Format(Fll .11 ,6X,F11 .11 ) . 

3.  Then,  it  reads  P cards  (P=3)  each  having  the  following 
parameters: 

a)  MTBF(J) , MTTR(J) , in  Format(F20.8,lX,F20.8). 

b)  A(J),  the  procurement  cost  for  component  J,  in 
Format(F10.2) . 

c)  B(J),  the  repair/replacement  cost  for  component  J,  in 
Format(F10.2) . 

d)  C(J),  the  cannibalization  cost  for  component  J,  in 
Format(F10.2) . 

4.  Then,  it  reads  the  revenue  values,  V(l),  and  V(2),  for 
having  1,  and  2 units  up,  respectively,  in 
Format(F10.2,lX,F10.2). 


After  this,  it  reads  the  sets  of  feasible  spare  allocations, 
M(1 ) , M ( 2 ) , and  M(3),  in  Format (II ,1X,I1 ,1X,I1 ). 

The  output  and  program  documentation  are  self-explanatory. 
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K T B F ( I ) = 0.80123000  MTTfU  1}=  0.1C000000 
MTBF(2)=  1.60200000  MTTR(2)=  0.12019000 
MTfcF ( 3 ) = 80.12800090  MTTR13)=  0.12019000 
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